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ABSTRACT 

The flowfield about a vertically-launched surface-to-air 
missile model at an angle of attack of 50 degrees and a 
Reynolds number of 1.1x10 s was investigated in a low-speed 
wind tunnel at the Naval Postgraduate School. The location 
and intensity of the asymmetric vortices in the wake of the 
missile model were determined and the vortices were displayed 
using planar velocity vector, total pressure coefficient, and 
vorticity plots. The model configuration tested was a body- 
only configuration (wings, strakes, and tails removed). Two 
flowfield conditions were treated: the nominal ambient wind 
tunnel condition and a grid-generated turbulence condition. 
Flow visualization was conducted and video-taped for both the 
body-only configuration and the winged configurations. The 
following conclusions were reached: 1) the addition of 
turbulence decreased the vorticity but did not significantly 
change the patterns of the plots; 2) the addition of 
turbulence reduced the vorticity more at eleven body diameters 
than at six body diameters; 3) compared to the body-only case, 
the vorticity is reduced for the "x" case but not for the "+" 
case for the turbulence condition; 4) flow visualization 
verified vortices movement away from the missile as the tested 
point was moved aft along the missile body. 
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I. 



INTRODUCTION 



A. BACKGROUND 

Missile development in recent years has emphasized the 
capability for operation at high angle of attack due to demand 
for higher maneuverability and special mission requirements. 
This trend is amplified by the increasing demand to minimize 
the overall wingspan. Several examples of missile design 
which involve high angles of attack are vertical launching, 
missile storage separation from aircraft, launch from an 
aircraft container, the one-man anti-tank weapon, and the air- 
to-air and air-to-ground short range missile. These are 
examples where high angles of attack can be reached at 
subsonic Mach numbers. The flight Mach number in some cases 
can even range from low subsonic to high supersonic. In such 
flight conditions it is important to understand the flowfield 
around the missile and the effects on aerodynamic 
characteristics. [Ref. 1] 

An example of such a missile is the U.S. Navy's 
vertically-launched surface-to-air missile (VLSAM) . Upon 
launch it enters the ocean's surface environment at low 
velocity where it may encounter significant crosswinds 
resulting in a high angle of attack situation. In addition, 
the VLSAM launch conditions may include some degree of 
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turbulence caused by interaction of the air movement above the 
ocean surface, the ship's motion, and superstructure 
interference. [Ref. 2] High angle-of-attack conditions may 
result in the formation of asymmetric vortices along the 
missile body. These vortices induce side forces on the 
missile and can affect its flight aerodynamics. A vertically- 
launched anti-submarine rocket was lost due to asymmetric 
vortices caused by its high angle of attack and nose asymmetry 
[Ref. 3]. 

This thesis is a part of the research conducted at the 
Naval Postgraduate School (NPS) to understand the effects of 
turbulence on a VLSAM model at high angle of attack. 
Initially, Roane developed the missile model and procedure to 
measure various flowfield turbulence effects [Ref. 4]. Rabang 
studied the asymmetric-vortex-induced side force by varying 
turbulence, angle of attack, and roll angle. He concluded 
that turbulence with length scales comparable to the nose- 
generated vortices increased the side force and decreased flow 
unsteadiness. He also verified that nose roll angle can cause 
significant changes in the side forces and that the addition 
of low-aspect wings did not significantly change the magnitude 
of the side forces. [Ref. 2] 

Lung mapped the location and shape of the asymmetric 
vortices, in a cross-plane six missile diameters from the 
nose, using a body-only configuration with and without 



2 



freestream turbulence, at 50 degree angle of attack. He 
concluded that the strength of the vortices weakened and 
became more diffuse with turbulence. [Ref. 5] 

Viniotis repeated the procedure at six diameters using two 
wing configurations: zero degree roll angle ("+" 
configuration) and 45 degree roll angle ("x" configuration). 

He concluded that the addition of turbulence decreased vortex 
strength and the addition of wings causes the vortices to move 
closer to the missile body. [Ref. 6] 

Johnson continued the study by mapping the flowfield in 
a cross-plane eleven diameters from the nose for the "+" and 
"x" configurations. He concluded that the addition of 
turbulence decreased the vortex strength with little change 
in the vortices' positions. He also noted that the vortex 
strength was greater for the "+" configuration than for the 
"x" configuration, but that the amount of change was much less 
with turbulence addition. [Ref. 7] 

The purpose of this research was to study the vortices for 
the body-only configuration in turbulent and non-turbulent 
flow at eleven missile diameters from the missile nose. In 
addition, flow visualization was conducted for all missile 
configurations. Comparisons with the results of Lung and 
Johnson will be made to correlate the vortices' 
characteristics at six diameters and eleven diameters for all 
missile configurations. 
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B. HIGH ANGLE OF ATTACK AERODYNAMICS 

An important phenomenon of high angle-of-attack 
aerodynamics is the flow separation on body, wing, and tail 
surfaces and the shedding of free vortices. The flow 
separates along a line on the body surface via an interaction 
of the external flow with the laminar or turbulent boundary 
layer of the body. The separation causes the formation of 
feeding sheets which transfer the vorticity generated in the 
boundary layer into the external flow, which then roll up on 
the lee side of the body to form vortices. [Refs. 1,8] 

A slender body of revolution experiences four distinct 
aerodynamic regimes as the angle of attack is increased from 
0 to 90 degrees. The transition from one regime to the next 
as a function of angle of attack is also dependent on nose 
shape, overall fineness ratio, crossflow Mach number, Reynolds 
number, roll angle, free stream turbulence, and surface 
roughness. [Refs. 9,10] 

Figure 1 shows schematics of the four flow patterns about 
an ogive-cylinder body and the corresponding characteristics 
of normal force and side force. At very low angles of attack 
(0 to 5 degrees) there is no noticeable boundary layer 
separation. In this regime the axial flow components dominate 
and the flow can be characterized as a classical potential 
flow field and an attached laminar or turbulent boundary 
layer. At intermediate angles of attack (5 to 20 degrees) 
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Figure 1. Flow Regimes [Ref. 2] 
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the boundary layer separates on the lee side of the body and 
rolls up into a symmetric vortex pair that originates from the 
apex of the body and extends to the base of the body. The 
vortices are steady with time and no side force or yawing 
moment is present. [Ref. 9] 

In addition to the primary symmetric vortex pair, a second 
pair of vortices exists with the same rotational direction as 
the primary pair, positioned near the separation line as shown 
in Figure 2 . The second interior vortices pair and the 
secondary vortices are induced by the primary vortices with 
rotation directions opposite to those of the primary vortices. 
The primary and the second vortices move away from the body 
and their cores become larger as the angle of attack 
increases. The second and secondary vortices are usually 
negligible in strength as compared to the primary vortices. 
[Ref . 11] 

In the next regime, 20 to 60 degrees, crossflow effects 
begin to dominate and the vortices may become asymmetric, 
producing a side force and yawing moment. The formation of 
asymmetric vortices is discussed in the next section. Lastly, 
at an angle of attack between 60 and 90 degrees, the crossflow 
completely dominates and the vortex shedding becomes 
unsteady, starting aft and moving forward as the angle of 
attack increases. The boundary layer is shed in the form of 
a Karman vortex street or random wake depending upon the 
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second vortex 




Figure 2. Cross-flow Streamlines [Ref. 11] 

Reynolds number, Mach number, body geometry, etc. [Refs. 
9,10,11] 

C. A8YMMETRIC VORTICES 

As the angle of attack transitions through the 20 to 60 
degree regime, the vortex flowfield on the lee side of slender 
bodies may become asymmetric, even without sideslip. These 
asymmetric vortices produce an undesirable side force and 
yawing moment which can make necessary a complex autopilot or 
degrade the mission effectiveness of the missile. The side 
forces may cause an unacceptable target miss distance, and the 
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yawing moments may reduce stability and controllability. 
[Refs. 9, 11] Reynolds number and Mach number affect the 
magnitude of the side force, but do not strongly influence the 
angle of attack for the onset of the side force. In general 
an increase in Mach number will decrease the magnitude of the 
side force. [Ref. 12,13] 

The cause of asymmetric vortices formation is not 
completely understood, but there are two generally accepted 
theories. The first theory is that boundary-layer-induced 
asymmetry at the flow separation points causes the vortex 
flowfield to become asymmetric. The maximum vortex-induced 
normalized side force will occur in the critical Reynolds 
number regime where it is possible to have supercritial 
separation on one side of the body and subcritical separation 
on the other. [Refs. 14,6] The second theory suggests the 
principle cause of vortex flow asymmetry to be a hydrodynamic 
(inviscid) instability in the initially symmetric vortex 
formation. These vortices, which increase in strength with 
angle of attack, interact with the surrounding potential 
flowfield to form the asymmetric characteristic. [Refs. 12,15] 

A vortex-switching phenomenon has been observed in which 
the vortex pattern rapidly switches from an almost symmetric 
to a highly asymmetric configuration [Ref. 13]. This 
switching phenomenon may result from a fluctuation in the 
separation characteristics on the missile, or may relate to 
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a second inviscid solution in the leeward flowfield [Ref. 16]. 

The large "steady" side forces observed by many investigators 
indicate that the detached vortices fluctuate about some mean 
asymmetric configuration [Ref. 13]. 

The missile fore-body has proven to be the dominating 
component concerning flow separation [Ref. 1]. The addition 
of wings does, however, cause the vortices to move closer to 
the body. This results in an increase of the effective angle 
of attack for the formation of the asymmetric vortices. [Ref. 
18] The addition of strakes to a missile body greatly 
increases the lift and generates additional vortices which may 
alter the effect of the asymmetric vortices. The addition of 
tails has little effect on the flowfield and side forces of 
the forebody. Rabang has shown that the missile model retains 
the side force and length scale effects (described in the next 
section) of the induced vortices with the addition of strakes 
and wings. [Ref. 2] 

D. TURBULENCE 

Turbulence represents the presence of random, short 
duration variations in a flowfield with a given mean velocity. 
When considering the effects of turbulence on a body in a 
flowfield, it is important to compare the scale of the body 
to that of the turbulence. [Ref. 2] 
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Turbulence intensity measures the relative magnitude of 
high frequency velocity fluctuations in the flowfield. 
Greater intensities correspond to more kinetic energy in the 
turbulent flow. The turbulence intensity, T u , is 
mathematically expressed by the ratio of the streamwise root- 
mean-square (rms) velocity fluctuation, u', to the mean 
velocity component of the flowfield, U». [Ref. 2] 

T u =u'/U°o (1) 

Turbulence length scales are a measure of the time- 
averaged size of the fluid disturbance eddies. The turbulence 
length scale to body size ratio may determine the amount to 
which the turbulence affects the missile in the flowfield. 
Comparing the length scale to the missile length or missile 
diameter gives three possible results. For a length scale 
much greater than the body, the effect on the missile is 
similar to that of a steady-state flowfield. For a length 
scale on the same order as the body, the flowfield is non- 
steady and may cause rolling, pitching and yawing motion of 
the body. [Ref. 2] For a length scale much smaller than the 
body, most noticeably smaller than the missile diameter, the 
small-scale turbulence may affect the boundary layer 
development and flow separation over the missile body. [Ref. 
19] 

Since the length scale is an average of the disturbances, 
turbulence length scales of varying magnitude can be found in 
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a flowfield. These magnitude differences can be explained by 
the cascade effect. The turbulence eddies in the flowfield 
experience a strain in one direction which affects the scale 
of motion in the other directions as required by the 
conservation of angular momentum. The strain causes the 
eddies to break up into smaller scale and decreasing energy 
disturbances. The "cascade" of energy of turbulent motion 
continues to smaller and smaller scales until viscosity 
dissipates the smallest eddies. As the turbulence decreases 
in size and energy, the individual intensities decrease at a 
faster rate. These changes account for the dominance of the 
larger scale turbulence. [Ref. 20] 
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II. EXPERIMENT AND PROCEDURES 



A. OVERVIEW 

Pressure measurements were recorded in the VLSAM model 
wake to study the behavior of the asymmetric vortices. The 
body-only missile model was positioned at 50 degrees angle of 
attack for freestream conditions with and without generated 
turbulence in the Naval Postgraduate School low-speed wind 
tunnel. The pressure measurements were converted to total and 
planar velocity components, total pressure coefficients, and 
vorticity which were plotted to observe the vortices' 
characteristics . 

Figure 3 shows the survey grid in an x-y plane which was 
perpendicular to the freestream flow 11 missile diameters aft 
from the missile nose, measured along the missile body. The 
survey grid measured 3.5 inches by 5.0 inches. Incremental 
step size was 0.25 inches, yielding 15 columns in the 
horizontal direction each having 21 points in the vertical 
direction. The grid was centered on the missile body and 
started 1.50 inches away from the missile body due to physical 
constraints of the sting mounting arm and probe. 
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Figure 3. The Planar Survey Grid [Ref. 5] 
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B . APPARATUS 



Primary equipment used were the NPS wind tunnel, 
turbulence-generating grid , VLSAM model, three-axis traverser, 
five-hole pressure probe, data acquisition system, and data 
reduction/display software. 

1. Wind Tunnel and Turbulence Grids 



The low-speed, single-return, horizontal-flow wind 
tunnel is located in Halligan Hall at NPS (Figure 4). 




Figure 4. Naval Postgraduate School Wind Tunnel [Ref. 21] 
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The tunnel is powered by a 100-horsepower electric 
motor with a three-blade variable pitch fan and a four-speed 
Dodge truck transmission. Stator blades aft of the fan help 
to straighten the flow. Turbulence is reduced by two wire 
mesh screens upstream of the settling chamber and corner 
turning vanes. The test section is 45 inches by 32 inches 
with corner lighting enclosures that reduce the section area 
from 10 ft 2 to 9.88 ft 2 . A reflection plane in the test 
section decreases the useable height to 28 inches. The 
tunnel's contraction ratio is approximately 10:1. The tunnel 
was designed for test section velocities up to 290 ft/sec, 
operating at atmospheric pressure. A remote-controlled, 
flush-mounted turntable in the test section allows changing 
of model angle of attack in a horizontal plane. Wind tunnel 
temperature was measured in the settling chamber with a dial 
thermometer. Dynamic pressure was determined by the pressure 
difference between the static pressure in the settling chamber 
and in the test section by a water micro-manometer. The 
settling chamber and test section each have four static taps 
connecting to the manometer via a common manifold. [Ref. 21] 
The measured pressure difference, measured in centimeters of 
water by the manometer, is converted to the test-section 
reference velocity by equation (2) . 

U M ={ (2*2.046*cm H 2 0) / (K*p ) } 1/2 (2) 
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Where : 



U M = measured velocity (ft/sec) 

2.046 = conversion factor 

cm H 2 0 = manometer reading in cm of H 2 0 
K = wind tunnel calibration factor (for specific 

grid) 

p = air density (lb^ft 2 ) 

The wind tunnel calibration factor, K, corrects for 
the fact that the actual dynamic pressure is slightly 
different than the measured static pressure difference between 
the static pressure in the settling chamber and in the test 
section. The value of the calibration factor was determined 
by plotting the actual dynamic pressure measured by a pitot- 
static tube in the test section versus the measured pressure 
difference. [Ref. 21] The relationship was found to be linear 
where the slope of the curve is the tunnel calibration factor. 
The value of K was determined to be 0.8891 for the no-grid 
configuration and 1.6545 for grid three. A large wall 
disturbance due to the grid frame accounts for the change in 
K. [Ref. 6] 

Four sguare-mesh, sguare-bar, biplanar grids were 
designed and built by Roane to generate turbulence of varying 
intensities and length scales. Each grid, three constructed 
of wood and one of wire, was mounted in a wooden frame shown 
in Figure 5. [Ref. 4] The grid, placed 73 inches upstream of 
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Figure 5. Turbulence-Generating Grids 
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the model pivot (Figure 6) , creates nearly isotropic 
homogeneous turbulence. Table 1 summarizes measured 
turbulence intensities and estimated length scales determined 
by Roane. 

2 . VLSAM Model and Mount 

The VLSAM model is a cruciform tail-control missile 
with low aspect ratio wings. The nose, wings, strakes, and 
tails are detachable from a hollow cylindrical body which can 
be rotated in 45 degree increments. Figure 7 portrays the 
VLSAM specifications. [Ref. 2] The model was rigidly mounted 
in the test section using a sting mount attached to a rotating 
arm which allows the angle of attack to be changed (Figure 6) . 

3. Traverser and Five-Hole Pressure Probe 

The Velmex 8300 3-D traverser used three 
microcomputer-controlled stepping motors mounted on top the 
test section to allow for manual or computer operated movement 
of an attached five-hole pressure probe (Figure 8) . Traverser 
movement was controlled using software (PPROBE) and a 
microcomputer . 

Pressure measurements at each data point were acquired 
using a nulling five-hole pressure probe attached to the 
traverser (Figure 9) . Pressure sample averages for each port 
were recorded at each data point and used to determine 
velocity, pitch angle, total pressure, total pressure 
coefficient, static pressure, and static pressure coefficient. 
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Figure 6. Planview of VLSAM Model with Pressure Probe and 
Grid in the Test Section of the Wind Tunnel (not drawn to 
scale) [Ref. 2] 



Ci id 


Intensity 

(percent) 


Length Scale 

(in.) 


Turbulence/ 
Model Dia. 


Dynamic Pressure 
(Ib/ft 1 ) 


One 


3.31 


1.84 


1.05 


15.35 


Two 


2.78 


1.56 


0.89 


14.88 


Three 


1.88 


1.08 


0.62 


16.38 


hour 


0.47 


0.27 


0.15 


15.61 


None 


0.23 


- 


- 


15.85 



Table 1. Grid Turbulence Parameters [Ref. 5] 
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Total length 
Base diameter 
Length/diameter ratio 
Ogive nose length 
Ogive/diameter ratio 
Wing span/root chord 
Tail span/root chord 
Center of pressure 



22.85 in. 

1.75 in. 

13.06 
4.0 in. 

2.29 

3.13 in./ 13.55 in. 

5.50 in./ 1.70 in. 

13.5 inches aft of nose tip (approx.) 




8.510 



NAVAL POSTGRADUATE SCHOOL 
SURFACE-TO-AIR 
MISSILE MODEL 




. > 30 * 

— 7 

^ i 



0.250 - 



3.500 - 



3.084 1 






13 . 454 - 



l.eeo 



1.750 



15 . 350 - 




- 1.700 



Figure 7. Specifications of VLSAM Model [Ref. 2] 



20 




Figure 8. 3-D Traversing Assembly 




Figure 9. Five-Hole Pressure Probe [Ref. 22] 
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Probe port PI measured stagnation pressure. Ports P2 and P3 
were static ports connected to a portable manometer. Flow yaw 
angle was read directly by rotating the probe until P2 and P3 
were equal. P4 and P5 were used with calibration curves to 
determine pitch angle. Velocities were found using the static 
and stagnation pressures as described previously. [Ref. 22] 

4. Data Acquisition 

A 48-port scanivalve attached to the five-hole probe 
allowed the pressures to be measured and converted to voltages 
by one transducer. A Hewlett-Packard data acquisition system 
consisting of a relay multiplexer, digital multimeter, relay 
actuator, and software enabled system control from a 
microcomputer. [Refs. 6,7,24] 

C. EXPERIMENTAL CONDITIONS 

Test conditions used by Roane [Ref. 4], Lung [Ref. 5], 
Viniotis [Ref. 6], and Johnson [Ref. 7] were repeated to allow 
for data comparison. 

(1) Turbulence grid three was used to generate freestream 
turbulence. Rabang found that grid three generated vortex- 
length-scale turbulence which had the largest side force 
effect on the missile. [Ref. 2] 

(2) A Reynolds number of 1.1*10 5 was obtained using 
reference pressure differences of 7.2 cm H 2 0 for no grid and 
10.0 cm H 2 0 for grid three. These settings were determined 
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by Roane [Ref. 3] to insure comparable test section velocities 
for each tunnel configuration. These corrections are 
necessary because the pressure difference measured by the 
manometer includes an error resulting from disturbances caused 
by the screen. The calibration constants determined by Roane 
were used to calculate the actual dynamic pressure. 

(3) Three body configurations were used. (Figure 10) 

- Body A: wings and tails at 0 degree roll angle 

- Body B: no wings or tails 

- Body C: wings and tails at 45 degree roll angle 



Body B was used for pressure measurements and all 
configurations were used for flow visualization. 
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Figure 10. VLSAM Model Body Configurations [Ref. 2] 



(4) The nose roll angle was set at position eight. 
Rabang [Ref. 2] found that the missile model's nose roll angle 
affected the side force coefficient, as suggested by previous 
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investigators. He determined that position eight produced the 
largest side force. 

(5) Missile angle of attack was set at 50 degrees, 
determined by Rabang to be the angle of attack for the maximum 
side force. [Ref. 2] 

(6) The data plane was at eleven body diameters from the 
nose corresponding to that used by Johnson. [Ref. 7] 

(7) Wind tunnel temperature was not allowed to increase 
by more than 20 degrees during a data run. When the 
temperature reached the limit the run was stopped until the 
temperature cooled, at which point the testing resumed. 

D . SOFTWARE 

PPROBE (Appendix A) is a program written by Kindelspire 
[Ref. 25] and modified by Lung [Ref. 5] to control traverser 
movement. The operator may choose manual or computer control, 
the size of plane to be measured, and the step size for each 
movement. 

CALP (Appendix B) is a program correlating transducer 
output voltage and pressure using a calibration manometer. 
The program was executed prior to and upon completion of the 
data collection run. CALP's output was converted to a slope- 
intercept equation using linear regression on a hand-held 
calculator. The two equations were averaged and the results 
used in CONVERT to calibrate the scanivalve transducer output. 
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CONVERT (Appendix C) is a program written by Lung [Ref. 
5] which reads the PPROBE data file and computes total 
velocity, x and y components of velocity, pitch angle, total 
pressure, total pressure coefficient, static pressure, and 
static pressure coefficient for each data point. Velocities, 
pitch angle, and pressures were calculated using calibration 
curves provided by the probe manufacturer. Pressure 
coefficients were calculated using ambient room temperature 
and tunnel calibration factor (K) , and were non- 
dimensionalized by the dynamic pressure. Equations (3) and 
(4) were used to calculate the static and total pressure 
coefficients, respectively. 



Ss = ( p s l - p ! >/ 0 


(3) 


C p<=( p ,r p t >/Q 


(4) 



Where: 

C ps = Static pressure coefficient 
C pt = Total pressure coefficient 
Q = Freestream dynamic pressure 
P s = Freestream static pressure 
P t = Freestream total pressure 
P sl = Local static pressure 
P tl = Local total pressure 

VORTIC is a program written by Johnson [Ref. 7] to compute 
the dimensionless vorticity using the velocity output from 
CONVERT using equation (5) . (Appendix D) 
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V = { d (v/U°o) / d (x/ d) } - { 8 (u/U°o) / d (y/ d) } 



(5) 



Where : 

V = vorticity 

v = y component of the velocity 
x = streamwise coordinate distance 
u = x component of the velocity 
y = transverse coordinate distance 

TECPLOT is a commercial graphics program used to generate 
the velocity vector, pressure coefficient contour, and 
vorticity contour plots. 

E. EXPERIMENTAL PROCEDURE 

The following experimental procedure was used to collect 
data . 

(1) The data acquisition equipment and wind tunnel were 
energized. Ambient tunnel temperature and pressure were 
recorded. The wind tunnel speed was monitored using the 
appropriate manual manometer value corresponding to the tunnel 
configuration previously discussed. 

(2) Pressure calibration was performed and the results 
recorded . 

(3) PPROBE was executed. The probe was positioned to the 
start point using manual control. Computer control was 
selected followed by grid size, increment distance and input 
file name. The probe was physically rotated until the 
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difference between P 2 and P 3 was approximately zero as 
indicated by the portable digital manometer. The yaw angle 
was read from the probe and entered on the computer. Each of 
the five pressure ports was automatically sampled ten times, 
averaged, and displayed on the computer screen. The values 
displayed were in units of volts and subsequently converted 
to pressure units. The P 1 measurement had to be positive or 
a very small negative (> -0.5 volts) value to ensure the probe 
was facing the correct direction. The difference between P 2 
and P 3 values had to be less than 0.1 volts, which was 
approximately a 6.4% error of the dynamic pressure for the no- 
grid configuration and an 8.6% error of the dynamic pressure 
for the grid #3 configuration (the difference between P 2 and 
P 3 would ideally be zero) . If either requirement was not 
satisfied, the probe was adjusted and the point resampled. 
At the completion of each column of data points PPROBE would 
automatically store the x-y position and corresponding average 
values and move the traverser to the next column. This 
procedure was repeated for the remaining data plane. 

(4) CALP program was executed, the results were recorded, 
and the calibration curve was calculated as previously 
described . 

(5) CONVERT and VORTIC programs were executed to obtain 
necessary data for plotting results. 
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(6) Output files from CONVERT and VORTIC were converted 
to input file format for TECPLOT. Velocity vectors, pressure 
coefficient, and vorticity contours were then plotted. 

F. FLOW VISUALIZATION EQUIPMENT 

The flow visualization system was designed for the NPS 
low-speed wind tunnel by Chlebanowski [Ref. 26] and modified 
by Sommers [Ref. 27]. The system is depicted in Figure 11. 

Smoke generated external to the tunnel by a portable 
fog/smoke machine was piped to a 10-gallon plenum chamber 
which served to dampen surges of smoke from the machine and 
allow the smoke to cool slightly. The smoke exited the plenum 
chamber and entered the contraction section of the tunnel via 
flexible hoses and copper tubing. The smoke then entered the 
freestream tunnel air three feet ahead of the missile model. 

A laser sheet was used to illuminate a thin plane on the 
leeward side of the missile, allowing the flow phenomenon to 
be observed and recorded. The laser source was a class four 
Spectra-Physics , five watt. Argon-ion laser. The laser 
emitted a 1.25 mm blue-green beam through a Newport 
Corporation F-LFI laser fiber to a fiber optic handset, beam 
expander, and piano-cylindrical lens. The handset, expander, 
and lens were mounted to a manual traverse mechanism above the 
test section. The traverse allowed the laser sheet to be 
moved to any position along the missile. 
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Figure 11. Flow Visualization System [Ref. 27] 

The flow visualization runs were video-taped using a 
color video camera, 3/4-inch VCR, and color television 
monitor. Still photographs were taken of selected tape frames 
by the NPS Educational Media Department directly from the 
monitor screen. 
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III. RESULTS AND DISCUSSION 



The following sections discuss the velocity vector plots, 
total pressure coefficient contour plots, vorticity contour 
plots, and flow visualization photographs for the missile 
model configuration B (body-only) at eleven diameters. The 
results produced without using a turbulence grid are labeled 
"OB" and with turbulence grid three are labeled "3B" . 

The plots generated from Lung's data are at a survey grid 
plane six diameters (10.5 inches) aft from the missile nose. 
Lung also used the body-only configuration, but used grid one 
for turbulence generation. 

The plots view the missile as it was mounted in the tunnel 
when observed in the direction of the tunnel freestream. All 
measurements are in inches (indicted by ") , followed by the 
equivalent missile model base diameters (1.75 inches and 
indicated by d) , and will be referenced from the leeward side 
of the missile body and the centerline. The terms "above" and 
"below" are relative to the missile centerline. "Away" 
indicates the distance from the center of the missile body. 



A. CONFIGURATION OB 

The velocity vector plot (Figure 12) indicates the 
location of three vortices and shows the location of a saddle 
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VELOCITY VECTORS 
(OB) 




Figure 12. Velocity Vectors-Conf iguration OB 
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point. The two larger vortices are asymmetric with the flow 
in opposite directions. The upper one is located about 1" 
(0.57d) above and 2 1/8" (1.2 Id) away. The lower outside one 
is located 1 1/2" (0.86d) below and 6 1/4" (3.57d) away. The 
lower inner one is located 7/16" (0.25d) below and 2 1/4" 
(1.29d) away. The saddle point is located 5/8" (0.36d) below 
and 3 1/8" (1.79d) away. 

The total pressure coefficient plot (Figure 13) shows one 
large center of pressure loss located on the centerline and 
4 1/4" (2.34d) away. The center does not coincide with any 
core from the velocity vector plot. The contours change from 
C pt level five to level two (-0.29 to -0.82) in 3/4" in the 
upper half and in 1" in the lower half. This larger gradient 
in the upper half indicates a greater vortex strength. 

The vorticity plot (Figure 14) shows three centers. The 
largest upper one is 5/16" (0.18d) above and 5 1/4" (3.00d) 
away. The smaller upper one, which is weaker, is 5/16" 
(0.18d) above and 4 1/2" (2.57d) away. The lower one is 1 
3/16" (0.79d) below and 4 7/8" (2.79d) away. The zero value 
alternates above and below the missile centerline with no 
apparent pattern. The upper vorticity cores have positive 
vorticity. The lower one has negative vorticity and is 
smaller in magnitude than either of the upper two. 

The vorticity plot from Lung (Figure 15) shows two 
centers. The upper one core is 3/4" (0..43d) above and 2 7/8" 
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TOTAL PRESS. COEFF. 




Figure 13. Total Pressure Coefficient-Configuration OB 
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Figure 14. Vorticity-Configuration OB 




Figure 15. Vorticity-Conf iguration OB (six diameters) 
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(1.64d) away. The lower one is 5/8" below (0.36(1^ and 3 1/4" 
(1.86d) away. The upper core has the largest magnitude and 
the greatest gradient. There also appear to be two centers 
at the inner edge of the survey grid, one above and one below 
the centerline. 

B. CONFIGURATION 3B 

The velocity vector plot (Figure 16) indicates the 
location of three vortices and a saddle point, similar to the 
OB vector plot. The upper one is located 7/8" (0.50d) above 
and 2 1/8" (1.2 Id) away. The larger lower one is located 1 
1/2" (0.86d) below and 6 3/8" (3.64d) away. The smaller lower 
one is located 5/8" (0.39d) below and 3 1/8" (1.79d) away. 
The saddle point is located 5/8" (0.36d) below and 3 1/8" 
(1.79d) away. The two larger vortices are asymmetric with 
opposite flow direction. 

The total pressure coefficient plot (Figure 17) shows one 
large center of loss located 5/16" (0.18d) above and 4 5/8" 
(2.64d) away. The center does not coincide with any core from 
the velocity vector plot. The contours change from C pt level 
six to level two (-0.28 to -0.86) in 3/4" (0.43d) in the upper 
half and in 1 1/4" (0.71d) in the lower half. This larger 
gradient change indicates a greater vortex strength in the 
upper half. 
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Velocity Vectors-Configuration 3B 
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Figure 17. Total Pressure Coefficient-Configuration 3B 
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The vorticity plot (Figure 18) shows two centers. The 
upper one is 1/2" (0.29d) above and 5 1/2" (3.14d) away. The 
lower core is 1 1/8" (0.64d) and 5 1/4" (3.00d) away. The zero 
value is roughly along the missile centerline. The upper core 
has positive vorticity. The lower one has negative vorticity 
and is smaller in magnitude than the upper core. 

The vorticity plot from Lung (Figure 19) represents 
turbulence added by use of grid #1. The plot shows one center 
1/2" (0.29d) above and 3" (1.71d) away. A second center is 
indicated in the lower half. The upper core has a much higher 
magnitude and gradient than the lower one. 

C. COMPARISON BETWEEN NO TURBULENCE AND TURBULENCE 

The addition of turbulence on the scale of the vortices 
did not affect the overall shapes of the plots. The effect 
on the velocity vector plot was to decrease the magnitudes, 
but the directions and core positions were approximately the 
same. The addition of turbulence did not change the 
magnitudes of the total pressure coefficient contours, but the 
gradient was much tighter and the center of the core moved up 
5/16" (0.18d). The vorticity strength was less and the 
gradients were greatly reduced with turbulence. The vorticity 
plots from Lung showed similar changes between the no 
turbulence and turbulence conditions at six diameters. 
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Figure 18. Vorticity-Conf iguration 3B 
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Figure 19. Vorticity-Conf iguration OB (six diameters) 
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These results indicate that the added turbulence smoothes 



and steadies the flow, and results in a decrease in vorticity. 
This decrease in vorticity for added turbulence corresponds 
to Rabang ' s findings that the side force decreases with added 
turbulence for the body-only configuration [Ref. 2]. 

D. COMPARISON BETWEEN BODY-ONLY CONFIGURATIONS 

Comparison of the vorticity plots for the body-only 
configuration at six diameters and eleven diameters for the 
no turbulence condition shows a higher vorticity strength at 
eleven diameters, but the gradients are approximately equal. 
The velocity vector plots (Figures 12,20) show that the 
positions of the vortex cores are farther away from the 
missile body at eleven diameters although the relative 
positioning is about the same. The total pressure coefficient 
plots (Figures 13,21) show two tight cores and a larger 
gradient at six diameters compared to one large core at eleven 
diameters. The second core is expected to be outside the 
survey grid. 

For the turbulence condition, however, the vortex strength 
is much higher and the gradient is larger at six diameters. 
The velocity vector and total pressure coefficient plots show 
the same relative characteristics as the no turbulence 
condition . 

These results indicate that the added turbulence reduces 
the vorticity at both six diameters and eleven diameters, but 
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Figure 21. Total Pressure Coefficient-Configuration OB (six 
diameters) 
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the magnitude of the vorticity reduction is greater at eleven 
diameters . 

E. COMPARISON BETWEEN BODY-ONLY AND WINGED CONFIGURATIONS 

Comparison of the body-only ("OB" and "3B") configuration 
to the "x" ("OC" and "3C") and "+" ("OA" and "3A") 

configurations [Ref. 7] at eleven diameters shows the body- 
only configuration to have a strong resemblance to the "x" 
configuration and little resemblance to the "+" configuration 
for both the no turbulence and turbulence conditions. 
Comparisons of the velocity vector plots between "OC" and "3C" 
(Figures 22,23) and "OB" and "3B" (Figures 12,16) show the 
same positioning and flow directions of the vorticity cores 
and the saddle points. Comparison of the total pressure 
coefficient plots (Figures 13,17,24,25) show comparable 
coefficient magnitudes, but the body-only configurations have 
a much larger gradient. Comparison of the vorticity plots 
with the "x" configuration shows a higher vorticity for the 
body-only configuration with no turbulence (Figures 14,26) and 
comparable vorticity magnitudes for the turbulence condition 
(Figures 18,27). Comparisons with the "+" configurations 

(Figures 28,29) show that the body-only configuration 

magnitude is greater with no turbulence but is less for the 
turbulence condition. 
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Figure 23. Velocity Vectors-Conf iguration 3C 
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Figure 24. Total Pressure Coefficient-Configuration OC 
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Figure 25. Total Pressure Coefficient-Configuration 3C 
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Figure 26. Vorticity-Conf iguration OC 
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Figure 27. 



Vorticity-Conf iguration 3C 
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Figure 28. Vorticity-Conf iguration OA 
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Figure 29. Vorticity-Configuration 3A 
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These results indicate that the addition of wings, 
strakes, and tails in the "x" configuration tends to reduce 
the magnitude of the vorticity with no turbulence. Again, as 
in the body-only configuration, the addition of turbulence of 
the nose-generated vortex scale appears to steady the flow and 
results in a decrease in vorticity for the "x" configuration. 
For the " + " configuration the turbulence has less of an effect 
on the vorticity magnitude, indicating a dependence on roll 
angle also. In addition, the amount of decrease in vorticity 
is larger for the body-only configuration. 

F. COMPARISON WITH FLOW VISUALIZATION 

Comparison of the body-only configuration velocity vector 
plots at six and eleven body diameters with flow visualization 
photographs (Figures 30,31) show that the velocity vector plot 
cores correspond to the asymmetric vortices formed on the lee 
side of the missile. Figure 32 is a sketch illustrating the 
vortex core movement away from the missile as the distance 
from the missile nose is increased from six to eleven body 
diameters . 
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Figure 30. 



Flow Visualization (six diameters) 
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Figure 31. 



Flow Visualization (eleven diameters) 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



The flowfield about a vertically-launched surface-to-air 
missile (VLSAM) model at high angle of attack was studied. 
The body-only configuration at eleven body diameters was 
tested and compared to the body-only configuration at six body 
diameters and with the winged configurations at eleven body 
diameters. The angle of attack was set at 50 degrees with a 
Reynolds number of 1.1*10 5 for all data runs. All tests were 
conducted with both an ambient wind tunnel condition (no grid) 
and a condition with grid-generated turbulence on the scale 
of the nose-generated asymmetric vortices. The following 
conclusions were reached: 

(1) Addition of turbulence for the body-only 
configuration tended to smooth and steady the flow. The 
turbulence reduced the vorticity level, but did not 
significantly change the positions of the vortices. 

(2) The addition of turbulence had a greater effect on 
vorticity reduction at eleven body diameters than at six body 
diameters for the body-only configuration, based on the 
vorticity plots (Figures 14,15,18,19). 

(3) The body-only plots at eleven body diameters are very 
similar to the "x" configuration at the same point, but 
different from the "+" configuration. Without turbulence, the 
vorticity levels for the body-only case are much larger than 
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for the "x" and the " + " configurations. With turbulence, 

however, the comparisons show that the vorticity levels were 

not only reduced but were approximately equal between the 

body-only and "x" configurations. However, compared to the 
"+" configuration, the vorticity levels were lower. 
Apparently the effect of added turbulence on vorticity is much 
greater for the body-only configuration than the winged 

configurations. The winged cases are also complicated by a 
dependence on roll angle. Rabang noted the same results on 
the side force magnitudes. 

(4) The addition of wings had a significant effect on the 
reduction of the vorticity level without added turbulence. 
With turbulence, the addition of wings had a negligible 
effect on the vorticity level. Again, Rabang noted the same 
results on the side force magnitude. Apparently, the effects 
of added turbulence outweighs the effects of adding wings. 

(5) The flow visualization experiment verified that the 
asymmetric vortices move away from the missile as distance is 
increased from the missile nose along the body axis, in a vary 
non-linear manner. The velocity vector plots tended to match 
the flow visualization results for the location of the vortex 
cores . 
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Recommendations for future research are: 



(1) Conduct pressure measurements outside the present 
survey grid to determine characteristics of secondary vortices 
that may not have been located near the missile or farther 
away from the missile. 

(2) Attempt to decrease the vortex asymmetry by modifying 
the present missile nose geometry. 
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APPENDIX A. PPROBE PROGRAM 



1 DEF SEG: CLEAR , &HFEOO-.GOTO 4 'Begin PCIB Program Shell 

2 GOTO 1000 ' User program 

3 GOTO 900 ' Error handling 

4 I=&HFE00 ' Copyright Hewlett-Packard 1984,1985 

5 PCIB. DIR$=ENVIRON$( "PCIB") 

6 I$=PCIB. DIR$+"\PCIBILC . BLD" 

7 BLOAD I$,I 

8 CALL I (PCIB. DIR$, 1%, J%) :PCIB.SEG=I% 

9 IF J%=0 THEN GOTO 13 

10 PRINT "Unable to load."; 

11 PRINT " (Error #";J%?")" 

12 END 

13 ' 

14 DEF SEG=PCIB.SEG:O.S=5:C.S=10:I.V=15 

15 I . C=2 0 : L . P=2 5 : LD . FILE=3 0 

16 GET . MEM=3 5 : L . S=4 0 : PANELS=4 5 : DEF . ERR=50 

17 PCIB . ERR$=STRING$ (64,32) : PCIB . NAME$=STRING$ ( 16 , 3 2 ) 

18 CALL DEF. ERR(PCIB. ERR, PCIB. ERR$, PCIB. NAME$, PCIB. GLBERR) ; 
PCIB. BASERR=255 

19 ON ERROR GOTO 3 

20 J=-l 

21 I$=PCIB. DIR$+"\PCIB. SYN" 

22 CALL O.S(I$) 

2 3 IF PCIB . ERRO0 THEN ERROR PCIB . BAS ERR 
24 1=0 

2 5 CALL 

I . V ( I , READ . REGISTER, READ . SELFID , DEFINE , INITIALIZE . SYSTEM) 

2 6 IF PCIB . ERRO0 THEN ERROR PCIB . BASERR 

2 7 CALL 

I . V ( I , ENABLE . SYSTEM , DISABLE . SYSTEM , INITIALIZE , POWER . ON ) 

28 IF PCIB . ERRoO THEN ERROR PCIB. BASERR 

29 CALL I. V( I, MEASURE, OUTPUT, START, HALT) 

30 IF PCIB. ERRO0 THEN ERROR PCIB. BASERR 

3 1 CAL 

I . V ( I , ENABLE . INT . TRIGGER , DISABLE . INT . TRIGGER , ENABLE . OUTPUT , 
IS ABLE. OUTPUT) 

3 2 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

33 CALL I. V( I, CHECK. DONE, GET. STATUS, SET. FUNCTION, SET. RANGE) 

34 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

35 CALL I. V( I, SET. MODE, WRITE. CAL, READ. CAL, STORE. CAL) 

3 6 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

37 CALL I. V(I, DELAY, SAVE. SYSTEM, J,J) 

38 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

39 1=1 

40 CALL I. V( I, SET. GATETIME, SET. SAMPLES, SET. SLOPE, SET. SOURCE) 
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41 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

42 CALL I . C ( I , FREQUENCY , AUTO • FREQ , PERIOD , AUTO . PER) 

43 IF PCIB . ERR<>0 THEN ERROR PCIB. BASERR 

44 CALL I. C(I, INTERVAL, RATIO, TOTALIZE, R100MILLI) 

4 5 IF PCIB . ERRoO THEN ERROR PCIB. BASERR 

46 CALL I . C ( I , R1 , RIO , R100 , R1KILO) 

47 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

48 CALL I . C ( I , R IOMEGA, RIO OMEGA , CHAN .A, CHAN . B) 

49 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

50 CALL I . C (I , POSITIVE , NEGATIVE , COMN, SEPARATE) 

51 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

52 1 = 2 

53 1 = 3 

54 CALL I .V(I, ZERO. OHMS, SET. SPEED, J,J) 

55 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

56 CALL I . C ( I , DCVOLTS , ACVOLTS , OHMS , R200MILLI ) 

57 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

58 CALL I . C ( I , R2 , R2 0 , R2 00 , R2KILO) 

59 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

60 CALL I . C ( I , R2 OKILO , R2 00 KILO , R2MEGA , R2 OMEGA) 

61 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

62 CALL I . C ( I , AUTOM ,R2 . 5,R12 . 5 , J) 

63 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

64 1=4 

6 5 CALL 

I . V ( I , SET . COMPLEMENT , SET . DRIVER , OUTPUT . NO . WAIT , ENABLE . HANDSH 
AKE) 

66 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

6 7 CALL 

I . V ( I , DISABLE . HANDSHAKE , SET . THRESHOLD , SET . START . BIT , SET . NUM . 
BITS) 

68 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

69 CALL I. V(I, SET. LOGIC. SENSE, J,J,J) 

7 0 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

71 CALL I. C( I, POSITIVE, NEGATIVE, TWOS, UNSIGNED) 

72 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

73 CALL I . C ( I , OC , TTL, R0 , R1 ) 

74 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

75 CALL I.C(I,R2,R3,R4,R5) 

76 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

77 CALL I . C ( I , R6 , R7 , R8 , R9 ) 

78 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

79 CALL I.C(I,R10,R11,R12,R13) 

8 0 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

81 CALL I.C(I,R14,R15,R16,J) 

82 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

83 1 = 6 

8 4 CALL 

I . V ( I , SET . FREQUENCY , SET . AMPLITUDE , SET . OFFSET , SET . SYMMETRY ) 

85 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

86 CALL I.V(I, SET. BURST. COUNT, J,J,J) 
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87 IF PCIB . ERRoO THEN ERROR PCIB. BASERR 

88 CALL I . C ( I , SINE , SQUARE , TRIANGLE , CONTINUOUS ) 

89 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

90 CALL I. C( I, GATED, BURST, J, J) 

91 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

92 1=7 

9 3 CALL 

I . V ( I , AUTOSCALE , CALIBRATE , SET . SENSITIVITY , SET . VERT . OFFSET) 

94 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

9 5 CALL 

I . V ( I , SET . COUPLING , SET . POLARITY , SET . SWEEPSPEED , SET . DELAY ) 

96 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

9 7 CALL 

I . V ( I , SET . TRIG . SOURCE , SET . TRIG . SLOPE , SET . TRIG . LEVEL , SET . TRIG 
.MODE) 

98 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

9 9 CALL 

I . V ( I , GET . SINGLE . WF , GET . TWO . WF , GET . VERT . INFO , GET . TIMEBASE . IN 
FO) 

100 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

10 1 CALL 

I . V ( I , GET . TRIG . INFO , CALC . WFVOLT , CALC . WFTIME , CALC . WF . STATS ) 
102 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

10 3 CALL 

I . V ( I , CALC . RISETIME , CALC . FALLTIME , CALC . PERIOD , CALC . FREQUENCY ) 
104 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

10 5 CALL 

I . V ( I , CALC . PLUSWI DTH , CALC . MINUSWIDTH , CALC . OVERSHOOT , CALC . PRE 
SHOOT) 

106 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

10 7 CALL 

I . V ( I , CALC . PK . TO . PK , SET . TIMEOUT , SCOPE . START , MEASURE . S INGLE . WF ) 

108 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

109 CALL I. V( I, MEASURE. TWO. WF, J, J, J) 

110 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

111 CALL I . C ( I , R10NANO , R100NANO , R1MICRO , R10MICRO) 

112 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

113 CALL I . C (I , R100MICRO, R1MILLI , R10MILLI , R100MILLI) 

114 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

115 CALL I . C ( I , R1 , R10 , R2 ONANO , R2 00NANO) 

116 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

117 CALL I . C ( I , R2MICRO , R2 OMICRO , R2 00MICRO , R2MILLI ) 

118 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

119 CALL I . C (I , R20MILLI , R200MILLI , R2 , R20) 

120 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

121 CALL I . C ( I , R5 ONANO , R50 ONANO , R5MICRO , R5 OMICRO) 

12 2 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

123 CALL I . C ( I , R50 OMICRO , R5MILLI , R50MILLI , R500MILLI ) 

124 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

125 CALL I . C ( I , R5 , R5 0 , CHAN . A , CHAN . B ) 

12 6 IF PCIB. ERROO THEN ERROR PCIB. BASERR 



63 



127 CALL I . C ( I , EXTERNAL, POSITIVE , NEGATIVE , AC) 

128 IF PCIB . ERRoO THEN ERROR PCIB . BASERR 

129 CALL I . C ( I , DC , TRIGGERED, AUTO . TRIG, AUTO . LEVEL) 

130 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

131 CALL I. C (I , XI, XI 0, STANDARD, AVERAGE) 

132 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

133 1=8 

134 CALL I. V( I, OPEN. CHANNEL, CLOSE. CHANNEL, J,J) 

135 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

136 CALL C.S 

137 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

138 I$=PCIB . DIR$+" \PCIB . PLD" 

139 CALL L . P ( I $ ) 

14 0 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

141 I$=" DMM . 01" : 1=3 :J=0: K=0 : L=1 

142 CALL DEFINE ( DMM .01,I$,I,J,K,L) 

14 3 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

144 I$="Func.Gen. 01": 1=6: J=0:K=1:L=1 

145 CALL DEFINE (Func . Gen .01,I$,I,J,K,L) 

14 6 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

147 I$="Scope. 01": 1=7 : J=0:K=2 :L=1 

148 CALL DEFINE (Scope .01,I$,I,J,K,L) 

14 9 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

150 I$="Counter . 01" : 1=1 : J=0 : K=3 : L=1 

151 CALL DEFINE ( Counter. 01,I$,I,J,K,L) 

152 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

153 I$="Dig. In. 01" : 1=4 : J=0 : K=4 : L=1 

154 CALL DEFINE ( Dig .In.01,I$,I,J,K,L) 

155 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

156 I$="Dig. Out . 01" : 1=4 : J=1 : K=4 : L=1 

157 CALL DEFINE ( Dig . Out .01,I$,I,J,K,L) 

158 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

159 I $= "Re lay . Act . 01" : 1=8 : J=0 : K=5 : L=1 

160 CALL DEFINE (Relay . Act .01,I$,I,J,K,L) 

161 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

162 I$="Relay . Mux . 01" : 1=2 : J=0 : K=6 : L=1 

163 CALL DEFINE (Relay . Mux. 01,I$,I,J,K,L) 

164 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

800 I $=ENVIR0N$ ( " PANELS " ) + " \PANELS . EXE " 

801 CALL L. S (1$) 

899 GOTO 2 

900 IF ERR=PCIB. BASERR THEN GOTO 903 

901 PRINT "BASIC error #";ERR;" occurred in line " ; ERL 

902 STOP 

903 TMPERR=PCIB. ERR: IF TMPERR=0 THEN TMPERR=PCIB. GLBERR 

904 PRINT "PC Instrument error #";TMPERR;" detected at line 
" ; ERL 

905 PRINT "Error: "?PCIB.ERR$ 

906 IF LEFT$ ( PCIB. NAME$ , 1) <>CHR$ ( 32 ) THEN PRINT "Instrument: 
" ; PCIB. NAME $ 

907 STOP 
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N 



908 COMMON PCIB. DIR$ , PCIB. SEG 

909 COMMON LD. FILE , GET . MEM, PANELS , DEF . ERR 

910 COMMO 
PCIB. BASERR, PCIB. ERR, PCIB. ERR$ , PCIB. NAME$ , PCIB. GLBERR 

911 COMMON 

READ. REGISTER, READ. SELFID, DEFINE , INITIALIZE . SYSTEM , ENABLE . SY 
STEM , DISABLE . SYSTEM , INITIALIZE , POWER . ON , MEASURE , OUTPUT , START 
, HALT , ENABLE . INT . TRIGGER , DISABLE . INT . TRIGGER , ENABLE . OUTPUT , D 
ISABLE . OUTPUT , CHECK . DONE , GET . STATUS 912 COMMON 
SET . FUNCTION , SET . RANGE , SET . MODE , WRITE . CAL, READ . CAL, STORE . CAL 
, DELAY , SAVE . SYSTEM , SET . GATETIME , SET . SAMPLES , SET . SLOPE , SET . SO 
URCE , ZERO . OHMS , SET . SPEED , SET . COMPLEMENT , SET . DRIVER , OUTPUT . NO 
.WAIT, ENABLE. HANDSHAKE, DISABLE. HANDSHAKE 913 COMMON 

SET . THRESHOLD, SET . START . BIT , SET . NUM. BITS , SET . LOGIC . SENSE , SET 
. FREQUENCY , SET . AMPLITUDE , SET . OFFSET , SET . SYMMETRY , SET . BURST . C 
OUNT , AUTOSCALE , CALIBRATE , SET . SENSITIVITY , SET . VERT . OFFSET , SET 
. COUPLING, SET. POLARITY, SET. SWEEPSPEED 914 COMMON 

SET . DELAY , SET . TRIG . SOURCE , SET . TRIG . SLOPE , SET . TRIG . LEVEL, SET . 
TRIG . MODE , GET . SINGLE . WF , GET . TWO . WF , GET . VERT . INFO , GET . TIMEBAS 
E . INFO , GET . TRIG . INFO , CALC . WFVOLT , CALC . WFTIME , CALC . WF . STATS , C 
ALC. RISETIME, CALC. FALLTIME, CALC. PERIOD 915 COMMON 

CALC. FREQUENCY, CALC. PLUSWIDTH, CALC. MINUSWIDTH, CALC. OVERSHOOT 
, CALC . PRESHOOT , CALC . PK. TO . PK , SET . TIMEOUT , SCOPE . START , MEASURE 
. SINGLE. WF, MEASURE. TWO. WF, OPEN. CHANNEL, CLOSE . CHANNEL 916 

COMMON 
FREQUENCY , AUTO . FREQ , PERIOD , AUTO . PER , INTERVAL , RATI O , TOTALI Z E , 
R100MILLI , R1 , RIO , R100 , R1KILO , R10MEGA, R100MEGA , CHAN . A, CHAN . B, 
POSITIVE , NEGATIVE , COMN , SEPARATE , DCVOLTS , ACVOLTS , OHMS , R2 OOMI L 
LI,R2,R20,R200, R2 KI LO , R2 OKI LO , R2 0 OKI LO 917 COMMON 

R2MEGA , R2 OMEGA , AUTOM , R2 . 5 , R1 2 . 5 , POSITIVE , NEGATIVE , TWOS , UNS IG 
NED, OC, TTL, RO , R1 , R2 , R3 , R4 , R5 , R6 , R7 , R8 , R9 , RIO , R1 1 , R12 , R13 , R14 
, R1 5 , R1 6 , SINE , SQUARE , TRIANGLE , CONTINUOUS , GATED , BURST , R1 ONANO 
, RIOONANO , R1MICRO , RIOMICRO , R1 OOMICRO 918 COMMON 

R1MILLI , R10MILLI , R100MILLI , R1 , RIO , R2 ONANO , R200NANO, R2MICRO, R 
2 OMICRO, R2 OOMICRO, R2MILLI , R20MILLI , R200MILLI , R2 , R20 , R50NANO, 
R50 ONANO , R5MICRO , R50MICRO , R5 OOMICRO , R5MILLI , R50MILLI , R500MIL 
LI, R5,R50, CHAN. A, CHAN. B, EXTERNAL, POSITIVE 919 COMMON 

NEGATIVE , AC , DC , TRI GGERED , AUTO .TRIG, AUTO . LE VE L , X 1 , X 1 0 , STAN DAR 
D , AVERAGE 

920 COMMON 

DMM . 01 , Func . Gen. 01 , Scope. 01 , Counter. 01 , Dig. In . 01 , Dig. Out . 01 , 
Relay . Act . 01 , Relay .Mux. 01 999 'End PCIB Program Shell 
1000 REM This step initialzes the HP system 
1010 CLS 

1020 OPTION BASE 1 

1 0 3 0 DIM 

P ( 5 ) , PA (50,5) , PP (50,5) ,XPT(50) , YPT(50) ,X(50),Y(50) , YAW (50) 
1040 REM 

1050 CALL INITIALIZE. SYSTEM (PGMSHEL. HPC) 

1060 REM 
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'RELAY MUX 



'RELAY 



107 0 REM SET FUNCTIONON THE ' DMM ' , 

ACTUATOR ' 

1080 REM 

1090 CALL SET. FUNCTION (DMM. 01, DCVOLTS) 

1100 CALL S ET. RANGE (DMM. 01, AUTOM) 

1110 CALL DISABLE. INT. TRIGGER (DMM. 01) 

1120 CALL ENABLE. OUTPUT (RE LAY. MUX. 01) 

1130 CALL ENABLE. OUTPUT (RE LAY. ACT. 01) 

1140 REM ************* PROGRAM TRAVERSE ******************* 
1150 REM 

1160 REM OPEN THE COM PORT AND INITIALIZE THE MOTOR 

SETTINGS 

1170 OPEN "coral : 1200 , n, 8 , 1 , rs , cs, ds , cd" AS #1 
1180 REM SET MOTOR DEFAULT VALUES 

1190 DATA 2000,2000,2000,2,2,2,0.000125,0.000125,0.000125 
1200 READ V1,V2,V3,R1,R2,R3,C1,C2,C3 

1210 REM DEFINE CHARACTERS FOR DATA REDUCTION ALGORITHM 
1220 RN2$="RENAME A : RAW . DAT " 

1230 HEAD1$ = " # X Y PI P2 P3 

P4 P5 YAW " 

1240 FORMAT$= "## ##.## ##.## ###.### ###.### ###.### 

###.### ###•### ###•##" 

1250 PRINT 

1260 PRINT "*******************************************" 

1270 PRINT "** USER MUST SELECT 'CAPS LOCK' FUNCTION **" 

1280 PRINT M ******************************************* n 
1290 REM DISPLAY MOTOR DEFAULT SETTINGS 



1300 PRINT " 

1310 PRINT " 
SETTINGS: " 

1320 PRINT " 

1330 PRINT " 
STEPS/SEC " 2 ) " 

1340 PRINT " 
INCHES" 

1350 PRINT " 

1360 PRINT 
1370 PRINT "NOTE!! 
POSITION BEFORE" 
1380 PRINT " 



♦★A********************************" 

INITIALIZED VALUES FOR ALL MOTOR 

VELOCITY = 1000 STEPS/SEC" 

RAMP (MOTOR ACCELERATION) = 2 (6000 

DEFAULT INCREMENTAL UNITS ARE 

******************************* t * * * * •' 

USE MANUAL CONTROL TO INITIALIZE PROBE 
SELECTING COMPUTER CONTROLLED MOVEMENT. 



1390 PRINT 

1400 INPUT "MANUAL CONTROL OR COMPUTER CONTROL (ENTER 'MAN' 
or ' CP ' ) " ; CON$ 

1410 IF CON$="CP" THEN 3490 

1420 REM OPTION TO CHANGE DEFAULT SETTINGS OF VELOCITY OR 
ACCELERATION RAMP 
1430 PRINT 
1440 PRINT 

1450 PRINT " DO YOU WANT TO CHANGE THE VELOCITY OR 
ACCELERATION RAMP" 
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DEFAULT SETTINGS? (Y or N) " 



1460 PRINT " 

1470 PRINT 

1480 PRINT "IF 'NO', THIS PROGRAM WILL THEN LET YOU DEFINE 
THE" 

1490 PRINT "DISTANCE YOU WANT TO MOVE (IN INCHES). IF 'YES'," 
1500 PRINT "YOU CAN CHANGE ANY OR ALL OF THE DEFAULT SETTINGS 
FOR ANY MOTOR." 

1510 PRINT 
1520 PRINT 
1530 PRINT 

154 0 INPUT "DO YOU WANT TO CHANGE ANY OF THE DEFAULT SETTINGS? 
(Y or N) " ; D$ 

1550 IF D$=" Y" THEN 1590 
1560 IF D$="N" THEN 2220 
1570 REM 

1580 REM **** OPERATOR SELECTED MOTOR VARIABLES ***** 

1590 PRINT 
1600 PRINT 

1610 INPUT "WHICH DEFAULT VALUE? (ENTER ' 1 ' FOR VELOC OR '2' 

FOR ACCEL RAMP) " ; L 

1620 ON L GOTO 1690,1930 

1630 PRINT "DO YOU WANT TO CHANGE THE DEFAULT VELOCITY? (Y OR 
N) " 

1640 INPUT V$ 

1650 IF V$="Y" THEN 1690 

1660 PRINT "DO YOU WANT TO CHANGE THE DEFAULT ACCELERATION 
RAMP? (Y or N) " 

1670 IF R$=" Y" THEN 1990 
1680 IF R$="N" THEN 1450 
1690 PRINT 
1700 PRINT 

1710 INPUT "WHICH MOTOR VELOCITY DO YOU WISH TO CHANGE? (1,2, 
or 3 ) " ; J 

1720 ON J GOTO 1730,1830,1880 
1730 PRINT 
1740 PRINT 

1750 INPUT "ENTER DESIRED VELOCITY OF MOTOR #1";V1 
1760 PRINT 
1770 PRINT 
1780 PRINT 

1790 PRINT "DO YOU WANT TO CHANGE VELOCITY OF ANOTHER MOTOR? 
(Y OR N) " 

1800 INPUT V$ 

1810 IF V$=" Y" THEN 1690 
1820 IF V$="N" THEN 1430 
1830 PRINT 
1840 PRINT 

1850 INPUT "ENTER DESIRED VELOCITY OF MOTOR 2";V2 
1860 PRINT 
1870 GOTO 1780 
1880 PRINT 
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1890 PRINT 

1900 INPUT "ENTER DESIRED VELOCITY OF MOTOR #3";V3 

1910 PRINT 

1920 GOTO 1780 

1930 PRINT 

1940 PRINT 

1950 INPUT "WHICH MOTOR ACCEL RAMP DO YOU WANT TO CHANGE? (1, 
2 , or 3 ) " ; K 

1960 ON K GOTO 1970,2060,2120 
1970 PRINT 
1980 PRINT 

1990 INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR #1";R1 
2000 PRINT 
2010 PRINT 

2020 PRINT "DO YOU WANT TO CHANGE THE ACCEL RAMP OF ANOTHER 
MOTOR? (Y or N)?" 

2030 INPUT RM$ 

2040 IF RM$=" Y" THEN 1930 
2050 IF RM$="N" THEN 1450 
2060 PRINT 
2070 PRINT 

2080 INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR #2";R2 

2090 PRINT 

2100 PRINT 

2110 GOTO 2000 

2120 PRINT 

2130 PRINT 

2140 INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR #3";R3 

2150 PRINT 

2160 PRINT 

2170 GOTO 2000 

2180 REM 

2190 REM DEFINE DISTANCE TO MOVE MOTOR 
2200 PRINT 
2210 PRINT 
2220 PRINT 

2230 REM INITIALIZE MOTOR INCREMENTS TO ZERO 

2240 11=0 

2250 12=0 

2260 13=0 

2270 PRINT 

2 2 8 OPR I NT" 

*************************************************************" 



2290 PRINT " ** 
* * » 

2300 PRINT " ** 
**" 

2310 PRINT " ** 
A MOTOR * * " 
2320 PRINT " ** 
MOTOR **" 



DEFINE WHICH MOTOR YOU WANT TO MOVE 



NOTE!!! A POSITIVE ('+') INCREMENT TO 
MOVES TRAVERSER AWAY FROM THAT PARTICULAR 
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2330 PRINT " ** 

**" 

2340 PRINT " ** — MOTOR #1 MOVES THE PROBE UPSTREAM AGAINST 

THE FLOW **" 

2350 PRINT " ** — MOTOR #2 MOVES THE PROBE TOWARD THE ACCESS 

WINDOW **" 

2360 PRINT " ** — MOTOR #3 MOVES THE PROBE VERTICALLY 

DOWNWARD **" 

2 3 7 OPR I NT" 

************************************************************* n 

2380 PRINT 
2390 PRINT 

2400 INPUT "WHICH MOTOR DO YOU WANT TO MOVE? (1,2, or 3)";L 
2410 ON L GOTO 2420,2680,2970 
2420 PRINT 
2430 PRINT 

2440 PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #1?" 

2450 PRINT " ********* (ENTER DISTANCE IN INCHES) *********" 
2460 INPUT II 
2470 PRINT 

2480 PRINT" *********************************" 

2490 PRINT 

2500 PRINT "SUMMARY OF OPERATOR INPUTS:" 

2510 PRINT " MOTOR #1 VELOCITY = ";V1 

2520 PRINT " ACCELERATION RAMP = ";R1 

2530 PRINT " INCREMENTAL DISTANCE = 

" ; II ; "INCHES" 

2540 PRINT" *********************************" 

2550 PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or 
N) " 

2560 PRINT 

2570 PRINT "ENTER 'N' TO START MOTOR MOVEMENT. ENTER 'Y' TO 
RETURN" 

2580 PRINT "TO VARIABLE SELECTION SUBROUTINE." 

2590 INPUT V$ 

2600 IF V$="Y" THEN 1430 
2610 GOSUB 3410 
2620 PRINT 

2630 PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y or N) ?" 
2640 INPUT C$ 

2650 IF C$=" Y" THEN 2220 
2660 IF C$="N" THEN 3260 
2670 PRINT 
2680 PRINT 

2690 PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #2?" 

2700 PRINT " ********* (ENTER DISTANCE IN INCHES) *********" 
2710 INPUT 12 
2720 PRINT 
2730 PRINT 

2740 REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 
2750 PRINT" *********************************" 
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2760 PRINT 

2770 PRINT "SUMMARY OF OPERATOR INPUTS:" 

2780 PRINT " MOTOR #2 VELOCITY = ";V2 

2790 PRINT " ACCELERATION RAMP = ";R2 

2800 PRINT " INCREMENTAL DISTANCE = 

12 INCHES" 

2810 PRINT" *********************************" 

2820 PRINT 
2830 PRINT 

2840 PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or 
N) " 

2850 PRINT 

28 60 PRINT "ENTER 'N' TO START MOTOR MOVEMENT. ENTER ' Y * TO 
RETURN" 

2870 PRINT "TO VARIABLE SELECTION SUBROUTINE." 

2880 INPUT V$ 

2890 IF V$="Y" THEN 1430 
2900 GOSUB 3410 
2910 PRINT 

2920 PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y or N) ?" 
2930 INPUT C$ 

2940 IF C$=" Y" THEN 2220 
2950 IF C$="N" THEN 3260 
2960 PRINT 
2970 PRINT 

2980 PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #3?" 

2990 PRINT " ********* (ENTER DISTANCE IN INCHES) *********" 
3000 INPUT 13 
3010 PRINT 
3020 PRINT 

3030 REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 
3040 PRINT" *********************************" 

3050 PRINT 

3060 PRINT "SUMMARY OF OPERATOR INPUTS:" 

3070 PRINT " MOTOR #3 VELOCITY = " ;V3 

3080 PRINT " ACCELERATION RAMP = " ;R3 

3090 PRINT " INCREMENTAL DISTANCE = 

" ?I3 ?" INCHES" 

3100 PRINT 

3110 PRINT" *********************************" 

3120 PRINT 
3130 PRINT 

3140 PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or 
N) " 

3150 PRINT 

3160 PRINT "ENTER 'N' TO START MOTOR MOVEMENT. ENTER ' Y ' TO 
RETURN" 

3170 PRINT "TO VARIABLE SELECTION SUBROUTINE." 

3180 INPUT V$ 

3190 IF V$=" Y" THEN 1430 
3200 GOSUB 3410 
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****************************11 

THE PROGRAM HAS ENDED . " 
****************************" 



3210 PRINT 
3220 PRINT 

3230 INPUT "DO YOU WANT TO INPUT ANOTHER MANUAL MOTOR MOVEMENT 
(Y or N) " ;M$ 

3240 IF M$=" Y" THEN 2210 
3250 PRINT 

3 260 PRINT "DO YOU WANT TO INPUT COMPUTER CONTROLLED MOTOR 
MOVEMENT?" 

3270 PRINT " ********* NOTE!!! ********* " 

3280 PRINT " ALL PREVIOUS MOTOR INCREMENT INPUTS HAVE BEEN 
ZEROIZED. " 

3290 PRINT "PROGAM WILL LET YOU CHOOSE MANUAL OR CP-CONTROLLED 
MOVEMENT . " 

3300 PRINT "***** (IF 'NO', THE PROGRAM WILL END). ****** 
3310 PRINT 

3320 INPUT "DO YOU WANT COMPUTER CONTROLLED MOTOR MOVEMENT (Y 
or N) " ;N$ 

3330 IF N$=" Y" THEN 3500 
3340 PRINT 
3350 PRINT 
3360 PRINT 
3370 PRINT " 

3380 PRINT " 

3390 PRINT " 

3400 END 

3410 REM ******* MOTOR MOVEMENT SUBROUTINE ********* 

3420 PRINT #1, »&» :PRINT #1, 

"E" ; "Cl=" ; Cl ; " :C2=" ;C2 ; " : C3=" ;C3 

3430 PRINT #1, "Il=" ; II ; " : Vl=" ; VI ; " :R1=" ;R1 ; 

3440 PRINT #1, " : I2=" ; 12 ; " : V2=" ; V2 ; " : R2=" ;R2 

3450 PRINT #1, "I3=" ; 13 ; " : V3=" ; V3 ; " : R3=" ;R3 ; " : 

3460 RETURN 

3470 REM ******************************************** 

3480 REM ******************************************** 

3490 PRINT 

3500 REM ******* COMPUTER CONTROLLED MOVEMENT ******* 

3510 PRINT 

3520 PRINT "THE PRESSURE DATA WILL BE WRITTEN TO FILES ON 
DRIVE ' A ' " 

3530 PRINT 

3540 PRINT "YOU WILL BE ASKED TO INPUT FILE NAMES FOR THESE." 
3550 PRINT 

3560 INPUT "IS A FORMATTED DISK IN DRIVE 'A'? PRESS 'ENTER' 
TO CONTINUE" ;D$ 

3570 PRINT 
3580 PRINT 
3590 PRINT 
3600 PRINT " 

3610 PRINT " 

3620 PRINT " 

3630 PRINT " 



***************************************** 
** note ! ! ! **" 

** COMPUTER CONTROLLED MOVEMENT **" 

** IS PROGRAMMED WITH A **" 
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** DEFAULTED NEGATIVE MOTOR INCREMENT **" 
** (i.e. MOTOR #3 WILL MOVE UPWARD **" 
** BY ENTERING A (+) DISTANCE) . **" 
***************************************** 



3640 PRINT " 

3650 PRINT " 

3660 PRINT " 

3670 PRINT " 

3680 PRINT 

3690 REM SET INITIAL 
POINTS TO ZERO 
3700 HT=0 
3710 WD=0 
3720 DIST=0 
3730 XPT=0 
3740 YPT=0 
3750 N=0 
3760 PRINT 
3770 PRINT 

3780 INPUT "WHAT IS THE DIMENSION 
YOU WANT TO MEASURE." ;WD,HT 
3790 PRINT 
38 00 INPUT "WHAT 
MOVE. " ; DIST 

3810 YPT=INT (HT /DIST) + 1 
3820 XPT=INT ( WD /DIST) + 1 
3830 N=XPT*YPT 
3840 PRINT 

3850 PRINT "THERE ARE ";XPT?" 
MEASURED " 

3860 PRINT 



MOVEMENT DISTANCE AND NUMBER OF DATA 



( X , Y ) (IN INCHES) THAT 



IS THE STEP (IN INCHES) THAT YOU WANT TO 



* " ; YPT ; " = ";N;" POINTS TO BE 



3870 

3880 

3890 

3900 

3910 

3920 

NAMES 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 



INPUT "ARE THE 
IF C$="N" THEN 
CLS 
N=XPT 

IF (N < 1) OR (N > 
REM *** GENERATING 



NUMBER 

3780 



OF POINTS IS OK. (Y OR N)";C$ 



99) GOTO 3780 
STRING STRING 



SEGMENTS FOR DATA FILE 



B$ = MID$ ( STR$ ( 1 ) , 2) : 
E$ = MID$ (STR$ (N) , 2) : 
X$ = "XXXXXX" 

EX$ = ".DAT" 

CLS 

"DATA FILES WILL 



REM 

REM 



** 
* * 



STRING 

ENDING 



NUMBER 

STRING 



ii •]_ ii 

NUMBER 



"N" 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
EXTENSION) " ;F2$ 
4040 PRINT 
4050 PRINT 
4060 IF LEN (F2$) 
4070 CLS 



BE INCREMENTED FROM:" 



(X$ + B$ + EX$) 



To 



(X$ + E$ + EX$) 



"ENTER DATA FILE NAME (6 CHARACTERS MAX — NO 



> 6 OR LEN (F2$) < 1 GOTO 4030 
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4080 PRINT N? "DATA FILES WILL BE GENERATED AND INCREMENTED 
AS FOLLOWS:" 

4090 PRINT 
4100 PRINT 

4110 PRINT ( F2 $ + B$ + EX$) ; " To (F2$ + E$ + EX$) 

4120 PRINT 
4130 PRINT 

4140 INPUT "ARE THE NUMBER OF POINTS AND FILE NAMES OK. (Y OR 
N)"; C$ 

4150 IF C$ = "N" GOTO 3780 

4160 IF C$ = "Y" GOTO 4180 

4170 GOTO 4140 

4180 CLS 

4190 PRINT 

4200 PRINT 

4210 REM SET INITIAL POSITION DATA 

4220 X ( 1) =-DIST 

4230 Y ( 1) =-DIST 

4240 FOR IX=2 TO XPT+1 

4250 X(IX)=0 

4260 NEXT IX 

4270 FOR JY=2 TO YPT+1 

4280 Y ( JY) =0 

4290 NEXT JY 

4300 FOR 1=1 TO XPT 

4302 11=0 

4304 12=0 

4306 13=0 

4310 FOR J=1 TO YPT 

4320 REM MOTOR CP-CONTROLLED MOTOR MOVEMENT 
4330 11=0 
4340 12=0 
4350 13=0 

4360 REM EACH POINT TAKE 10 TIMES READINGS 
4370 X ( 1+1) =X ( I ) +DIST 
4380 XPT(J) =X(I+1) 

4390 Y ( J+l) =Y ( J) +DIST 
4400 YPT (J) =Y (J+l) 

4405 INPUT " ADJUST THE WHEEL TO MAKE THE P2 =P3 , INPUT THE YAW 
ANGLE" ; YAW ( J) 

4408 PRINT 

4410 INPUT " PRESS 'ENTER' TO START THE MEASUREMENT" ; MOVE $ 
4420 REM 

4430 REM READ FIVE CHANNELS AND DISPLAY THE DATA 

4440 REM 

4450 STEPPER=4 

4460 SWITCH = 3 

4470 HOMER=8 

4480 DELAY 1 = .1 

4490 DELAY 2 = 1 

4500 REM SET THE S.V PORT TO #4 
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4510 FOR IL=1 TO 3 
4520 THYME = TIMER 

4530 CALL OUTPUT (RELAY . ACT . 01 , STEPPER) 

4540 CHKTIME = TIMER 

4550 IF CHKTIME < (THYME + DELAY 1 ) GOTO 4540 
4560 CALL OPEN . CHANNEL (RELAY . ACT . 01 , SWITCH) 

4570 CLS 
4580 NEXT IL 
4590 PRINT 

4600 PRINT " NOW IS POINT ";J 

4610 REM START MEASURE FROM PORT 4 TO PORT 8 
4620 FOR JJ=1 TO 5 

4630 CALL OUTPUT (RELAY . ACT . 01 , STEPPER) 

4640 CHKTIME = TIMER 

4650 IF CHKTIME < (THYME + DELAY2 ) GOTO 4640 
4660 REM EACH PORT SAMPLE 10 TIMES 
4670 FOR 11=1 TO 10 
4680 ROUT=l 

4690 CALL OUTPUT (RELAY . MUX . 0 1 , ROUT) 

4700 CALL MEASURE (DMM. 01, VOLTS) 

4710 PA (II , JJ ) =VOLTS 
4720 NEXT II 

4730 CALL OPEN. CHANNEL (RELAY .ACT . 01 , SWITCH) 

4740 IF JJ=5 THEN 4760 
4750 NEXT JJ 

4760 REM HOME THE S.V PORT TO #48 
4770 CALL OUTPUT (RELAY . ACT . 01 , HOMER) 

4780 CALL OPEN . CHANNEL (RELAY . ACT . 01 , HOMER) 

4790 REM 

4800 REM DISPLAY THE SAMPLE DATA 
4810 REM 

4820 PRINT HEAD1$ 

4830 FOR IS= 1 TO 10 

4840 PRINT USING 
FORMAT$ ; IS , XPT ( J) ,YPT(J) , PA ( IS , 1 ) , PA ( IS , 2 ) , PA ( IS , 3 ) , PA ( IS , 4 ) 

, PA ( IS , 5 ) , YAW ( J) 

4850 NEXT IS 
4860 REM 

4870 REM AVERAGE THE DATA 
4880 REM 

4890 FOR JA = 1 TO 5 

4900 TOTAL = 0 

4910 FOR IA = 1 TO 10 

4920 TOTAL = TOTAL + PA(IA,JA) 

4930 NEXT IA 

4940 AVERAGE = TOTAL /10 
4950 P ( JA) =AVERAGE 
4960 NEXT JA 
4970 PRINT 

4980 PRINT "THE AVERAGES ARE: " 

5000 PRINT HEAD1$ 
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5010 FOR JD=1 TO 5 
5020 PP(J, JD)=P(JD) 

5030 NEXT JD 

5040 PRINT USING 
FORMAT$ ; J , XPT ( J) ,YPT(J) , PP ( J , 1 ) , PP ( J , 2 ) , PP ( J , 3 ) , PP ( J , 4 ) , PP ( J 
,5) , YAW ( J) 

5045 PRINT 

5048 PRINT USING "THE NULLING ERROR IS +#.####" ;PP (J, 3 ) - 
PP ( J / 2 ) 

5049 PRINT 

5050 PRINT "DO YOU WANT RE-MEASURE AGAIN (Y / N) " 

5060 PRINT 

5062 PRINT "IF 'Y' WILL RE-SAMPLE AGAIN." 

5064 PRINT 

5070 INPUT "IF 'N' WILL MOVE THE TRAVERSER STEP UPWARD (WAIT 
7 SEC )";C$ 

5075 PRINT 

5080 IF C$=" Y" THEN 4405 
5082 IF C$="N" THEN 5090 
5084 GO TO 5070 
5090 IF J=YPT THEN 5160 
5100 REM 

5110 REM MOVE THE TRAVERSER STEP UPWARD. 

5120 REM 
5130 I3=-DIST 
5140 GOSUB 3410 
5150 NEXT J 

5160 REM*** STORE DATA BEFORE NEXT SAMPLE*** 

5170 OPEN "A:\RAW.DAT" FOR OUTPUT AS #2 
5180 PRINT #2 , HEAD1$ 

5190 FOR ID=1 TO YPT 

5200 PRINT #2 , USING 

FORMATS; ID, XPT (ID) , YPT (ID) , PP ( ID, 1 ) , PP ( ID , 2 ) , PP ( ID , 3 ) , PP ( ID, 
4 ) , PP ( ID , 5 ) , YAW (ID) 5210 NEXT ID 
5220 CLOSE #2 

5230 REM *** GENERATING INCREMENTED DATA FILE NAME 
5240 IF (I > 10) OR (I = 10) THEN 1$ = MID$ ( STR$ ( I ) , 2) 

5250 IF (I < 10) THEN 1$ = (MID$ ( STR$ ( 0 ) , 2) + MID$ ( STR$ ( I ) , 
2)) 

5260 FI2 $ = ( F2 $ + 1$ + EX$) 

5270 PRINT 

5280 PRINT " WRITING DATA FILE »; FI2$ 

5290 DF2$=RN2$+FI2$ 

5300 REM ** RENAME DATA FILE 
5310 SHELL DF2$ 

5320 REM 

5330 REM MOVE THE TRAVERSER TO THE NEXT SAMPLE POSITION 
5340 REM 
5350 PRINT 

5360 IF I=XPT THEN 5430 
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537 0 INPUT "THEN PRESS 'ENTER' FOR NEXT COLUMN SAMPLE ( 90 SEC) 
" ;MOVE$ 

5390 I2=-DIST 
5400 I 3=HT 
5410 GOSUB 3410 
5420 NEXT I 
5430 CLS 

5440 PRINT "ALL MOVEMENTS COMPLETE" 

5450 PRINT 
5460 PRINT 

5470 PRINT "YOU WANT TO REPOSITION TRAVERSER FOR ANOTHER 
MOVEMENT (Y OR N)?" 

5480 PRINT 

5490 PRINT "IF 'Y', THE PROGRAM WILL TAKE YOU TO MANUAL 
CONTROL SUBROUTINE." 

5500 PRINT "IF 'N', THE PROGRAM WILL END." 

5510 PRINT 

5520 INPUT "ANOTHER MOVEMENT" ;R$ 

5530 IF R$ = "Y" THEN 1370 
5540 IF R$ = "N" THEN 3370 
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APPENDIX B. CALP PROGRAM 



1 DEF SEG: CLEAR , &HFEOO:GOTO 4 'Begin PCIB Program Shell 

2 GOTO 1000 ' User program 

3 GOTO 900 1 Error handling 

4 I=&HFE00 ' Copyright Hewlett-Packard 1984,1985 

5 PCIB. DIR$=ENVIRON$( "PCIB") 

6 I$=PCIB. DIR$+"\PCIBILC. BLD" 

7 BLOAD I$,I 

8 CALL I (PCIB. DIR$, 1%, J%) :PCIB.SEG=I% 

9 IF J%=0 THEN GOTO 13 

10 PRINT "Unable to load."; 

11 PRINT " (Error #";J%;")" 

12 END 

13 ' 

14 DEF SEG=PCIB.SEG:O.S=5:C.S=10:I.V=15 

15 I . C=2 0 : L . P=2 5 : LD . FI LE=3 0 

1 6 GET . MEM=3 5 ; L . S=4 0 : PANELS=4 5 ; DEF . ERR=5 0 

17 PCIB.ERR$=STRING$ (64 , 32) : PCIB . NAME$=STRING$ ( 16 , 3 2 ) 

18 CALL DEF. ERR ( PCIB. ERR, PCIB . ERR$ , PCIB. NAME$, PCIB. GLBERR) : 
PCIB . BASERR=2 55 19 ON ERROR GOTO 3 

20 J=-l 

21 I$=PCIB. DIR$+"\PCIB. SYN" 

22 CALL O.S(I$) 

2 3 IF PCIB . ERRO0 THEN ERROR PCIB . BASERR 
24 1=0 

2 5 CALL 

I . V ( I , READ . REGISTER, READ. SELFID, DEFINE , INITIALIZE . SYSTEM) 

2 6 IF PCIB.ERRO0 THEN ERROR PCIB. BASERR 

2 7 CALL 

I . V ( I , ENABLE . SYSTEM , DISABLE . SYSTEM , INITIALI ZE , POWER . ON ) 

28 IF PCIB . ERRoO THEN ERROR PCIB. BASERR 

29 CALL I. V( I, MEASURE, OUTPUT, START, HALT) 

30 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

3 1 CAL 

I . V ( I , ENABLE . INT . TRIGGER, DISABLE . INT . TRIGGER, ENABLE . OUTPUT , 

I SABLE. OUTPUT) 32 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

33 CALL I. V( I, CHECK. DONE, GET. STATUS, SET. FUNCTION, SET. RANGE) 

34 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

35 CALL I. V( I, SET. MODE, WRITE. CAL, READ. CAL, STORE. CAL) 

3 6 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

37 CALL I.V(I, DELAY, SAVE. SYSTEM, J,J) 

38 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

39 1 = 1 

40 CALL I. V( I, SET. GATETIME, SET. SAMPLES, SET. SLOPE, SET. SOURCE) 

41 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

42 CALL I. C( I, FREQUENCY, AUTO. FREQ, PERIOD, AUTO. PER) 

4 3 IF PCIB. ERROO THEN ERROR PCIB. BASERR 
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44 CALL I. C(I, INTERVAL, RATIO, TOTALIZE, R100MILLI) 

4 5 IF PCIB.ERRO0 THEN ERROR PCIB . BASERR 

46 CALL I . C (I , R1 , RIO , R100 , R1KILO) 

47 IF PCIB . ERRO0 THEN ERROR PCIB. BASERR 

48 CALL I . C (I , R IOMEGA , RIO OMEGA, CHAN . A , CHAN . B) 

49 IF PCIB . ERRoO THEN ERROR PCIB. BASERR 

50 CALL I. C( I, POSITIVE,NEGATIVE,COMN, SEPARATE) 

51 IF PCIB. ERRO0 THEN ERROR PCIB. BASERR 

52 1=2 

53 1=3 

54 CALL I. V(I, ZERO. OHMS, SET. SPEED, J,J) 

55 IF PCIB. ERRO0 THEN ERROR PCIB. BASERR 

56 CALL I . C (I , DCVOLTS , ACVOLTS , OHMS , R200MILLI ) 

57 IF PCIB. ERRO0 THEN ERROR PCIB. BASERR 

58 CALL I . C (I , R2 , R20 , R200 , R2KILO) 

59 IF PCIB. ERRO0 THEN ERROR PCIB. BASERR 

60 CALL I . C ( I , R2 OKILO , R200KILO , R2MEGA , R2 OMEGA) 

61 IF PCIB. ERRO0 THEN ERROR PCIB. BASERR 

62 CALL I . C ( I , AUTOM , R2 . 5 , R1 2 . 5 , J ) 

63 IF PCIB. ERRO0 THEN ERROR PCIB. BASERR 

64 1=4 

6 5 CALL 

I . V ( I , SET . COMPLEMENT , SET . DRIVER , OUTPUT . NO . WAIT , ENABLE . HANDSH 
AKE) 66 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

6 7 CALL 

I . V ( I , DISABLE . HANDSHAKE , SET . THRESHOLD , SET . START . BIT , SET . NUM . 
BITS) 68 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

69 CALL I. V(I, SET. LOGIC. SENSE, J,J,J) 

70 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

71 CALL I. C( I, POSITIVE, NEGATIVE, TWOS, UNSIGNED) 

72 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

73 CALL I . C ( I , OC , TTL , R0 , R1 ) 

74 IF PCIB. ERRoO THEN ERROR PCIB. BASERR 

75 CALL I . C ( I , R2 , R3 , R4 , R5 ) 

76 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

77 CALL I . C ( I , R6 , R7 , R8 , R9 ) 

78 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

79 CALL I . C (I , RIO ,R11,R12,R13) 

80 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

81 CALL I.C(I,R14,R15,R16,J) 

82 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

83 1=6 

8 4 CALL 

I . V ( I , SET . FREQUENCY , SET . AMPLITUDE , SET . OFFSET , SET . SYMMETRY) 

85 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

86 CALL I. V(I, SET. BURST. COUNT, J,J,J) 

87 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

88 CALL I. C( I, SINE, SQUARE, TRIANGLE, CONTINUOUS) 

89 IF PCIB. ERROO THEN ERROR PCIB. BASERR 

90 CALL I. C(I, GATED, BURST, J,J) 

91 IF PCIB. ERROO THEN ERROR PCIB. BASERR 
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92 1=7 

9 3 CAL 

I . V ( I , AUTOSCALE , CALIBRATE , SET . SENSITIVITY , SET . VERT . OFFSET) 

94 IF PCIB.ERROO THEN ERROR PCIB.BASERR 
9 5 CALL 

I . V ( I , SET . COUPLING , SET . POLARITY , SET . SWEEPS PEED , SET . DELAY ) 

96 IF PCIB.ERROO THEN’ ERROR PCIB.BASERR 

9 7 CALL 

I . V ( I , SET . TRIG . SOURCE , SET . TRIG . SLOPE , SET . TRIG . LEVEL, SET . TRIG 
.MODE) 98 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

9 9 CALL 

I . V ( I , GET . SINGLE . WF , GET . TWO . WF , GET . VERT . INFO , GET . TIMEBASE . IN 
FO) 100 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

10 1 CALL 

I . V ( I , GET . TRIG . INFO , CALC . WFVOLT , CALC . WFTIME , CALC . WF . STATS ) 
102 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

10 3 CALL 

I . V ( I , CALC . RISETIME , CALC . FALLTIME , CALC . PERIOD , CALC . FREQUENCY) 
104 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

10 5 CALL 

I . V ( I , CALC . PLUSWIDTH , CALC . MINUSWIDTH , CALC . OVERSHOOT , CALC . PRE 
SHOOT) 106 IF PCIB.ERROO THEN ERROR PCIB.BASERR 
10 7 CALL 

I . V ( I , CALC . PK . TO . PK , SET . TIMEOUT , SCOPE . START , MEASURE . SINGLE . WF) 

108 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

109 CALL I . V ( I , MEASURE . TWO . WF , J , J , J ) 

110 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

111 CALL I . C ( I , R10NANO , R100NANO , R1MICRO , R10MICRO) 

112 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

113 CALL I . C ( I , R100MICRO , R1MILLI , R10MILLI , R100MILLI ) 

114 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

115 CALL I . C ( I , R1 , R10 , R2 ONANO , R200NANO) 

116 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

117 CALL I . C ( I , R2MICR0 , R2 OMICRO, R200MICRO , R2MILLI ) 

118 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

119 CALL I . C ( I , R20MILLI , R2 00MILLI , R2 , R2 0 ) 

12 0 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

121 CALL I . C (I , R50NANO, R500NANO, R5MICRO, R50MICRO) 

122 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

123 CALL I . C ( I , R50 OMICRO , R5MILLI , R50MILLI , R500MILLI ) 

124 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

125 CALL I.C(I,R5,R50, CHAN . A , CHAN . B ) 

12 6 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

127 CALL I. C( I, EXTERNAL, POSITIVE, NEGATIVE, AC) 

128 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

129 CALL I. C( I, DC, TRIGGERED, AUTO. TRIG, AUTO. LEVEL) 

130 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

131 CALL I. C( I, XI, XI 0, STANDARD, AVERAGE) 

132 IF PCIB.ERROO THEN ERROR PCIB.BASERR 

133 1=8 

134 CALL I. V( I, OPEN. CHANNEL, CLOSE. CHANNEL, J,J) 
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13 5 IF PCIB.ERROO THEN ERROR PCIB . BASERR 

136 CALL C.S 

137 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

138 I$=PCIB . DIR$+"\PCIB . PLD" 

139 CALL L. P ( 1$ ) 

14 0 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

141 I$="DMM. 01" : 1=3 : J=0 : K=0 : L=1 

142 CALL DEFINE ( DMM. 01,I$,I,J,K / L) 

14 3 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

144 I$="Func . Gen .01" : 1=6 : J=0 : K=1 : L=1 

145 CALL DEFINE ( FUNC . GEN .01,I$,I,J,K,L) 

14 6 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

147 I$="Scope . 01" : 1=7 : J=0 : K=2 : L=1 

148 CALL DEFINE ( SCOPE .01,I$,I,J,K,L) 

14 9 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

150 I$="Counter . 01" : 1=1 : J=0 : K=3 : L=1 

151 CALL DEFINE ( COUNTER .01,I$,I,J,K,L) 

152 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

153 I$="Dig. In. 01" : 1=4 : J=0:K=4 :L=1 

154 CALL DEFINE ( DIG . IN .01,I$,I,J,K,L) 

155 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

156 I$="Dig. Out .01" : 1=4 : J=1 : K=4 : L=1 

157 CALL DEFINE ( DIG . OUT .01,I$,I,J,K,L) 

158 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

159 I$="Relay . Act . 01" : 1=8 : J=0 : K=5 : L=1 

160 CALL DEFINE (RELAY . ACT .01 f I$ f I f J,K,L) 

161 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

162 I$="Relay .Mux. 01" : 1=2 : J=0 : K=6 : L=1 

163 CALL DEFINE (RELAY. MUX. 01,1$, I ,J,K,L) 

164 IF PCIB.ERROO THEN ERROR PCIB. BASERR 

800 I$=ENVIR0N$ ("PANELS" )+ "\PANELS.EXE" 

801 CALL L. S ( 1$) 

899 GOTO 2 

900 IF ERR=PC IB. BASERR THEN GOTO 903 

901 PRINT "BASIC error #";ERR;" occurred in line ";ERL 

902 STOP 

903 TMPERR=PCIB . ERR : IF TMPERR=0 THEN TMPERR=PCIB . GLBERR 

904 PRINT "PC Instrument error #";TMPERR;" detected at line 
" ; ERL 

905 PRINT "Error: ";PCIB.ERR$ 

906 IF LEFT$ ( PCIB . NAME$ , 1) <>CHR$ ( 32 ) THEN PRINT "Instrument: 
" ; PCIB . NAME$ 907 STOP 

908 COMMON PCIB . DIR$ , PCIB . SEG 

909 COMMON LD. FILE , GET . MEM, PANELS , DEF . ERR 

910 COMMON 
PCIB . BASERR , PCIB . ERR , PCIB . ERR$ , PCIB . NAME$ , PCIB . GLBERR 

911 COMMON 
READ . REGISTER, READ. SELFID, DEFINE , INITIALIZE . SYSTEM , ENABLE . SY 
STEM , DISABLE . SYSTEM , INITIALIZE , POWER . ON , MEASURE , OUTPUT , START 
, HALT , ENABLE . INT . TRIGGER, DISABLE . INT . TRIGGER , ENABLE . OUTPUT , D 
ISABLE . OUTPUT , CHECK . DONE , GET . STATUS 912 COMMON 
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SET . FUNCTION , SET . RANGE , SET . MODE , WRITE . CAL , READ . CAL , STORE . CAL 
, DELAY , SAVE . SYSTEM , SET . GATETIME , SET . SAMPLES , SET . SLOPE , SET . SO 
URCE , ZERO . OHMS , SET . SPEED , SET . COMPLEMENT , SET . DRIVER , OUTPUT . NO 
.WAIT, ENABLE. HANDSHAKE, DISABLE. HANDSHAKE 913 COMMON 

SET . THRESHOLD , SET . START . BIT , SET . NUM . BITS , SET . LOGIC . SENSE , SET 
. FREQUENCY , SET . AMPLITUDE , SET . OFFSET , SET . SYMMETRY , SET . BURST . C 
OUNT , AUTOSCALE , CALIBRATE , SET . SENSITIVITY , SET . VERT . OFFSET , SET 
.COUPLING, SET. POLARITY, SET. SWEEPSPEED 914 COMMON 

SET . DELAY , SET . TRIG . SOURCE , SET . TRIG . SLOPE , SET . TRIG . LEVEL, SET . 
TRIG . MODE , GET . SINGLE . WF , GET . TWO . WF , GET . VERT . INFO , GET . TIMEBAS 
E . INFO , GET . TRIG . INFO , CALC . WFVOLT , CALC . WFTIME , CALC . WF . STATS , C 
ALC. RISETIME, CALC. FALLTIME, CALC. PERIOD 915 COMMON 

CALC . FREQUENCY , CALC . PLUSWIDTH , CALC . MINUSWIDTH , CALC . OVERSHOOT 
, CALC . PRESHOOT , CALC . PK . TO . PK , SET . TIMEOUT , SCOPE . START , MEASURE 
. SINGLE . WF , MEASURE . TWO . WF , OPEN . CHANNEL , CLOSE . CHANNEL 916 

COMMON 
FREQUENCY , AUTO . FREQ , PERIOD , AUTO . PER , INTERVAL , RATIO , TOTALI Z E , 
R100MILLI , R1 , RIO , R100 , R1KILO, R1 OMEGA, RIO OMEGA, CHAN . A, CHAN . B, 
POSITIVE , NEGATIVE , COMN , SEPARATE , DCVOLTS , ACVOLTS , OHMS , R2 OOMIL 
LI,R2 ,R20,R200, R2 KI LO , R2 OKI LO , R2 0 OKI LO 917 COMMON 

R2 MEGA, R2 OMEGA, AUTOM, R2 . 5 , R12 . 5 , POSITIVE , NEGATIVE , TWOS , UNSIG 
NED, OC , TTL, RO , R1 , R2 , R3 , R4 , R5 , R6 , R7 , R8 , R9 , RIO , Rll , R12 , R13 , R14 
, R15 , R16 , SINE , SQUARE, TRIANGLE , CONTINUOUS , GATED, BURST, RIONANO 
, RIOONANO, R1MICRO, RIOMICRO, RIOOMICRO 918 COMMON 

R1MILLI , R10MILLI , R100MILLI , R1 , RIO , R20NANO , R2 OONANO , R2MICRO , R 
2 OMICRO , R2 OOMICRO , R2MILLI , R2 OMILLI , R2 OOMILLI , R2 , R2 0 , R50NANO , 
R500NANO, R5MICRO, R50MICRO, R500MICRO, R5MILLI , R50MILLI , R500MIL 
LI, R5,R50, CHAN. A, CHAN. B, EXTERNAL, POSITIVE 919 COMMON 

NEGATIVE , AC , DC , TRIGGERED , AUTO . TRIG , AUTO . LEVEL, XI , XI 0 , STANDAR 
D, AVERAGE 920 COMMON 
DMM . 01 , FUNC . GEN . 01 , SCOPE . 01 , COUNTER. 01 , DIG . IN . 01 , DIG . OUT . 01 , 
RELAY. ACT. 01, RELAY. MUX. 01 999 'End PCIB Program Shell 
1000 REM 

1010 REM This step initialzes the HP system 
1020 CLS 

1030 OPTION BASE 1 

1040 DIM P(10) ,PA(50,6) ,PP(50,6) ,XPT(40) ,CAL(40) 

1050 CALL INITIALIZE. SYSTEM ( PGMSHEL. HPC) 

1060 REM 

1070 REM All PC devices now have an initial state 
1080 REM Set function on the DMM and Relay MUX 
1090 REM 

1100 CALL SET. FUNCTION ( DMM . 01, DCVOLTS) 

1110 CALL SET. RANGE (DMM. 01, AUTOM) 

1120 CALL DISABLE. INT. TRIGGER (DMM. 01) 

1130 CALL ENABLE . OUTPUT ( RELAY. MUX. 01) 

1140 FORMAT$="## ##.#### ##.#### ##.#### ##.#### ##.#### 

##.####" 1200 FOR 1=1 TO 10 

1210 CAL ( I ) =0 . 0 

1220 NEXT I 

1510 REM 
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152 0 REM READ THE VOLTAGE OF 4 8TH CHANNEL AND DISPLAY THE 
DATA 

1530 REM 

1540 PRINT " CHOOSE 6 POINTS" 

1550 PRINT 

1550 PRINT "THE CALIBRATION WILL BE STORES IN 'CAL. DAT'" 

1560 REM 

1570 REM Begin sampling loop 
1580 REM 

1600 FOR J=1 TO 1 
1610 PRINT 

1630 FOR JJ=1 TO 6 

1631 INPUT "INPUT THE CALIBRATION PRESSURE" ; CAL (JJ ) 

1632 INPUT "PRESS 'ENTER' TO START MEASUREMENT" ;MOVE$ 

1640 FOR 11=1 TO 10 

1650 ROUT=l 

1660 CALL OUTPUT (RELAY. MUX. 01, ROUT) 

1670 CALL MEASURE (DMM. 01, VOLTS) 

1680 PA ( I I , J J ) =VOLTS 

1690 NEXT II 

1700 IF JJ=6 THEN 1740 

1730 NEXT JJ 

1740 REM 

1750 REM DISPLAY THE SAMPLE DATA 
1760 REM 

1780 FOR IS= 1 TO 10 

1790 PRINT USING 
FORMAT$ ; IS , PA ( IS , 1 ) , PA ( IS , 2 ) , PA ( IS , 3 ) , PA ( IS , 4 ) , PA ( IS , 5) , PA ( I 
S , 6 ) 1800 NEXT IS 

1810 REM 

1820 REM AVERAGE THE DATA 
1830 REM 

1840 FOR JA = 1 TO 6 

1850 TOTAL = 0 

1860 FOR IA = 1 TO 10 

1870 TOTAL = TOTAL + PA(IA,JA) 

1880 NEXT IA 

1890 AVERAGE = TOTAL /10 
1900 P ( J A ) =AVERAGE 
1920 NEXT JA 
1930 PRINT 

1940 PRINT "THE AVERAGE ARE: " 

2000 FOR JD=1 TO 6 
2010 PP(J, JD)=P(JD) 

2020 NEXT JD 

2055 PRINT USING 
FORMAT$ ; J , PP ( J , 1 ) ,PP(J,2) ,PP(J,3) ,PP(J,4) ,PP(J,5) ,PP(J,6) 2070 
PRINT 

2080 INPUT "DO YOU WANT RE-MEASURE AGAIN ? (Y / N)";C$ 

2090 IF C$=" Y" THEN 1580 

2101 REM*** STORE DATA BEFORE NEXT SAMPLE*** 
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2102 OPEN "A:\CAL.DAT" FOR OUTPUT AS #2 

2106 FOR ID=1 TO 6 

2107 PRINT #2 , USING FORMATS ; ID, PP (J , ID) , CAL ( ID) 

2108 NEXT ID 

2109 CLOSE #2 
2210 NEXT J 
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APPENDIX C. CONVERT PROGRAM 



************************************************************ 

***** 

* THIS PROGRAM CONVERTS THE VOLTAGE OF TRANSDUCER INTO 
PHYSICAL * 

* PRESSURE, VELOCITY, YAW ANGLE AND PITCH ANGLE. THOSE DATA ARE 
* 

* USED FOR PLOT PROGRAM LATER. 

* 

************************************************************ 

***** 

CHARACTER* 12 FNAME 
CHARACTER* 12 NAME 
CHARACTER* 12 FONAME 
CHARACTER* 2 A (50) 

CHARACTER*80 ST 

REAL K, INTR 

INTEGER COLS , RWS , DTPTS 





DATA A/ ' 01 ' , 


' 02 ' , 


' 03 ' 


1 

/ 


' 04 ' , 


, ' 05 ' , ' 06 ' , 


* 07 1 


1 

t 


'08 ' 


, '09' , 


* 


' 10 ' , 


'll' , 


' 12 ' 


1 

t 


' 13 ' , 


, ' 14 ' , ' 15' , 


' 16 1 


\ 

/ 


' 17 ' 


, '18', 


* 


' 19 ' , 


'20' , 


' 21 ' 


1 

/ 


' 22 ' , 


, '23 ' , '24 ' , 


' 25 1 


1 

/ 


'26' 


,' 21 ' , 


* 


' 28 ' , 


' 29 ' , 


' 30 ' 


1 

/ 


'31', 


, '32', '33', 


'34 1 


1 

/ 


' 35 1 


' , '36' , 


* 


' 37 ' , 


' 38 ' , 


' 39 ' 


/ 


'40' , 


, '41' , '42' , 


'43 ' 


1 

! 


'44 1 


1 , '45' , 


* 


'46' , 


'47 * , 


'48 ' 


1 

/ 


' 49 ' , 


, '50'/ 











WRITE 


(*, 


' ( A\) ' 


) ' # 


OF COLS 


(AWAY 


FROM 


MSL) = ' 


READ 


(*/ ' 


(15) ') 


COLS 










WRITE 


(*/ 


' ( A\) ' 


) ' # 


OF DATA 


PTS IN 


A COL (UP/ DOWN) = 


READ 


(*, ' 


(15) ') 


RWS 










WRITE 


(*, 


' ( A\) ' 


) ' DATA FILE 


NAME? 


(IE 


R001A2XX.DAT) 


READ 


(*, ' 


(A12) ' 


) NAME 










WRITE 


(*, 


' ( A\) ' 


) ' PI 


(F4.2) 


_ i 






READ 


(*, ' 


(F4.2) 


' ) PI 










WRITE 


(*, 


' ( A\) ' 


) ' PF 


(F4.2) 


— i 






READ 


(*, ' 


(F4.2) 


' ) PF 










WRITE 


(*, 


' ( A\) ' 


) ' TI 


(F3.1) 


_ i 






READ 


(*, ' 


(F3.1) 


') TI 










WRITE 


(*, 


' ( A\) ' 


) ' TF 


(F3.1) 


— i 






READ 


(*/ ' 


(F3.1) 


' ) TF 










WRITE 


(*, 


' ( A\) ' 


) ' K 


(F6.4) = 


i 






READ 


(*, ' 


(F6.4) 


') K 










WRITE 


(*/ 


' ( A\) ' 


) ' SLOPE FOR 


DELTAP 


( F9 . 


6) = ' 


READ 


(*, ' 


(F9 . 6) 


' ) SLOPE 








WRITE 


(*, 


' ( A\) ' ) ' INTERCEPT 


FOR DELTAP 


(F9.6) = ' 


READ 


(*, ’ 


(F9.6) 


' ) INTR 
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WRITE ( * , ' ( A\) 1 ) ' QM1 FACTOR (F4.2) = ' 

READ ( * , ' (F4.2) ') QM1FAC 
WRITE ( * , ' ( A\) ' ) ' X OFFSET = ' 

READ (*, ' (F5.2) ' ) XOFF 

WRITE ( * , ' ( A\) ' ) ' Y OFFSET = ' 

READ (*, • (F5.2) ' ) YOFF 

WRITE ( * , ' ( A\) ' >) ' OUTPUT FILE NAME = ' 

READ ( * , 1 ( A12 ) ' ) FONAME 

* CONVERT THE PRESSURE UNIT FROM inHg TO psf 

PATM= ( PI+PF ) *35.3631 
R= 1716. 5 
E=0 .0123 

T= (TI+TF) /2 . +4 60 
RO=PATM/ (R*T) 

DTPTS=RWS*COLS 

* OPEN A NEW FILE TO STORE THE REDUCED DATA 

OPEN ( 2 , FILE=FONAME , STATUS= ' NEW ' ) 

WRITE (2, 222) DTPTS 
222 FORMAT (15) 

* OPEN A SEQUENTIAL OF DATA FILE 

DO 20 1=1, COLS 
NAME (7:8) =A ( I ) 

FNAME=NAME 

OPEN ( 1 , FILE=FNAME) 

READ ( 1 , 100 , END=20) ST 
100 FORMAT (A6 5) 

15 READ ( 1 , 1000 , END=30 ) N0,X,Y,V1,V2,V3,V4,V5 ; BETA 

1000 FORMAT ( 12 , F7 . 2 , F6 . 2 , 5F9 . 3 , F8 . 2 ) 

* CONVERT THE VOLTAGE TO PRESSURE IN LBF/FT**2 

P1=DELTAP(V1, SLOPE, INTR) *2 . 0475+PATM 
P2=DELTAP (V2, SLOPE, INTR) *2 . 0475+PATM 
P3=DELTAP(V3 , SLOPE, INTR) *2 . 0475+PATM 
P4=DELTAP(V4 , SLOPE, INTR) *2 . 0475+PATM 
P5=DELTAP(V5, SLOPE, INTR) *2 . 0475+PATM 

* CALCULATE THE PITCH ANGLE IN DEGREES 

P= ( P4 -P5 ) / ( P1-P2 ) 

IF (P.GT.0.80) P=0 . 8 0 
IF (P.LT.-0.80) P=-0 . 8 0 
ALPHA=FPITCH (P) 

* CALCULATE THE VELOCITY IN FT/SEC 

YSLOP=FYSLOP (ALPHA) 

VELM=SQRT ( ( 2*YSLOP* ( P1-P2 ) )/ (RO*K) ) 

VEL=VELM* ( 1+E) 

* CALCULATE THE LOCAL DYNAMIC PRESSURE 

QM1=QM1FAC*2 . 0475/K 
QM=RO*VEL**2/2 . 

Q 1— QM1 * (1+2*E) 

Q=QM* ( 1+2*E) 

* CALCULATE THE YAW ANGLE IN DEGREES 

YAW=FYAW ( BETA+5 . 0 ) 

* CALCULATE THE VELOCITY COMPONENTS 
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BETAR= YAW*. 017453 

ALPHAR= (ALPHA-17. 942) *.017453 

VELY=VEL*SIN (ALPHAR) 

VELX=VEL*COS (ALPHAR) *SIN (BETAR) 

* CALCULATE THE TOTAL PRESSURE IN LBF/IN**2 

PTC=FPT (ALPHA) 

PT1=P1-Q*PTC 
PT=PT1/144 . 

CPT= (PT1-PATM-Q1) /Q1 

* CALCULATE THE STATIC PRESSURE IN LBF/IN**2 

PS1=PT1-Q 

PS=PS 1/144 

CPS= ( PS1-PATM) /Q1 

* WRITE VALUES TO OUTPUT FILE 

WRITE (2,2000) -X+XOFF , Y+YOFF , VEL, VELX , VELY , YAW , 

C ALPHA-17. 94 2, PT,CPT, PS, CPS 

2000 FORMAT (11F10. 3) 

GO TO 15 
30 CLOSE (1) 

20 CONTINUE 

CLOSE ( 2 ) 

STOP 

END 

***************************************************************** 

* THIS FUNCTION CONVERTS THE VOLTAGE TO PHYSICAL PRESSURE 

FUNCTION DELTAP (X, SLOPE, INTR) 

REAL INTR 

DELTAP=X*SLOPE+INTR 

END 

************************************************************* 

* THIS FUNCTION CALCULATES THE PITCH ANGLE 

FUNCTION FPITCH (X) 

FPITCH=3 . 759+53 .7568 *X-1. 308 5*X**2-1. 6583 *X**3 
* -0.8061*X**4+16.5115*X**5 

END 

* THIS FUNCTION CALCULATES THE VELOCITY PRESSURE COEFFICENT 

FUNCTION FYSLOP (X) 

IF ( X . LT . - 10 ) THEN 

FYSLOP=0 . 981-0 . 0102 *X- 3 . 000E-4*X**2-2 . 500E-6*X**3 
ELSE IF ( (X.GE.-10) .AND. (X.LE.10) ) THEN 
FYSLOP=0 . 98-0 . 006*X+2 . 000E-4*X**2 
ELSE 

FYSLOP=0. 9801-0. 0035*X-1. 143E-4*X**2+5. 833E-6*X**3 
END IF 
END 

************************************************************* 

* THIS FUNCTION CALCULATES THE YAW ANGLE 

FUNCTION FYAW(X) 

IF ( (X. GE . 0) .AND. (X . LE . 180) ) THEN 
FYAW=-X 



86 



ELSE 

FYAW=3 60-X 
END IF 
END 

************************************************************* 
* THIS FUNCTION CALCULATES THE TOTAL PRESSURE COEFFICIENT 
FUNCTION FPT(X) 

IF (X . LE .-30) THEN 
FPT=-0 . 01 

ELSE IF( (X.GT.-30) .AND. (X.LT.-20) ) THEN 
FPT=0 . 02+1 . 00E-3*X 

ELSE IF( (X.GE.-20) .AND. (X.LE.30) ) THEN 
FPT=0 
ELSE 

FPT=0 . 03-1 . 00E-3*X 
END IF 
END 
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APPENDIX D. VORTIC PROGRAM 



CHARACTER* 12 FNAME , OFNAME 

REAL VEL, VX (13,23) ,VY(13,23) ,VOR(13,23) ,VORX(13,23) 
REAL VELX, VELY , X , Y , DH , VORY (13,23) ,XC(13,23) ,YC(13,23) 
WRITE (*,' (A\) ' ) ' DATA FILE NAME ? ' 

READ ( * , ' ( A12 ) ' ) FNAME 

WRITE ( * , ' ( A\) ' ) ' OUTPUT FILE NAME ? ' 

READ ( * , ' ( A12 ) 1 ) OFNAME 

OPEN ( 3 , FILE=OFNAME , STATUS= ' NEW ' ) 

OPEN ( 2 , FILE=FNAME) 

10 READ ( 2 , 100 , END=20) X, Y, VEL, VELX, VELY 
100 FORMAT (5F10.3) 

C COMPUTE INDICES FOR ARRAYS 
J=INT ( (Y/.25J+12.0) 

I=INT ( (-X/.25) -1.0) 

C COMPUTE NON-DIMEN. VELOCITIES 
VX ( I , J ) =VELX/ VEL 
VY ( I , J ) =VELY/ VEL 
XC(I, J)=X 
YC ( I , J) =Y 
GO TO 10 

C COMPUTE NON-DIMEN. STEP SIZE 

C DELTA H=2 *GRID STEP DISTANCE/MISSILE DIAMETER 
20 DH= (2.0*.25)/1.75 

C COMPUTE VALUES FOR THE Y VORTICITY IN X DIR. ARRAY 
DO 3 0 J=1 , 2 3 
DO 40 1=1,13 

C IF STATEMENTS DEFINE BOUNDARIES 
IF(I.EQ.l) THEN 

VORY ( I , J ) = ( - 3 . * VY ( I , J ) +4 . * VY ( 1+ 1 , J ) -VY ( 1+ 2 , J ) ) / DH 
ELSE IF ( I . EQ . 15 ) THEN 

VORY ( I , J ) = ( 3 . * VY ( I , J ) -4 . *VY (1-1 , J) +VY ( 1-2 , J) )/DH 
ELSE 

VORY ( I , J) = (VY ( 1+1 , J) -VY ( 1-1 , J) ) /DH 
ENDIF 

40 CONTINUE 

30 CONTINUE 

C COMPUTE VALUES FOR THE X VORTICITY IN Y DIR. ARRAY 
DO 35 J=1 , 2 3 
DO 45 1=1,13 

C IF STATEMENTS DEFINE BOUNDARIES 
IF ( J . EQ . 1 ) THEN 

VORX (I,J)=(— 3. *VX ( I , J) +4 . *VX (I , J+l) -VX ( I , J+ 2 ) ) / DH 
ELSE IF ( J . EQ .21) THEN 

VORX ( I , J) = ( 3 . *VX ( I , J) -4 . *VX ( I , J-l ) +VX ( I , J-2 ) ) /DH 
ELSE 

VORX (I, J)=(VX(I, J+l) -VX(I, J-l) ) /DH 
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ENDIF 

45 CONTINUE 
35 CONTINUE 

C COMPUTE THE VORTICITY FOR EACH POINT 
DO 50 1=1,13 
DO 55 J=1 , 23 

VOR (I , J) =VORY ( I , J) -VORX ( I , J) 

WRITE (3,200) XC ( I , J) , YC(I,J), VOR(I,J) 
200 FORMAT (3F10.3) 

55 CONTINUE 
50 CONTINUE 
CLOSE (2) 

CLOSE (3) 

STOP 

END 
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